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ABSTRACT 

In the companion Paper XIX of this series we derive accurate total mass-to-light ratios 
{M /i) JAM = [M/L] {r — Rc) within a sphere of radius r — Rc centred on the galaxy, as well 
as stellar (M/L)stars (with the dark matter removed) for the volume-limited and nearly mass 
selected (stellar mass M^, ^ 6 x 10^ Mq) ATLAS'^'^ sample of 260 early-type galaxies (ETGs, 
ellipticals Es and lenticulars SOs). Here we use those parameters to study the two orthogonal 
projections (Mjam, <Je) and (Mjam , -R"'"'') of the thin Virial Plane (VP) {Mjam, <Je, n'r"") 
which describes the distribution of the galaxy population, where A/jam = L x (A//L)jam ~ 
The distribution of galaxy properties on both projections of the VP is characterized by (i) a 
boundary in the galaxy distribution, described by two power-laws, joined by a break at a char- 
acteristic mass A/jam ~ 3 x IO^^A/q, which corresponds to the minimum R^ and maximum 
stellar density, and (ii) a characteristic mass Mjam ~ 2 x lO^^Af© which separates a popula- 
tion dominated by fast rotator with disks at lower masses, from one dominated by quite round 
slow rotators at larger masses. The distribution of ETGs properties on the two projections of 
the VP tends to be constant along lines of constant Ce, or constant -y/ GM^,/ (SiR'^'^^) ~ ae 
respectively, and forms a continuous and parallel sequence with the distribution of spiral galax- 
ies. This applies to the dynamical (A//iy)jAM and to other indicators of the (A//L)pop of the 
stellar population, like H/3 and colour, as well as to galaxy concentration, which we show 
is tracing the bulge mass. A similar variation along contours of cTg is also observed for the 
mass normalization of the stellar Initial Mass Function (IMF), which was recently shown to 
vary systematically within the ETGs galaxy population. Our preferred relation has the form 
logio[(A//L),tars/(A^/i)saip] =a + bx logig (rje /130km s" ^ ) with a = -0.11 ± 0.01 and 
b = 0.36 ± 0.06. This trend implies a transition of the mean IMF from Kroupa to Salpeter in 
the interval log]^Q((Te/km s^^) ss 1.9 — 2.4 (or de w 80 — 260 km s~^ ) , with a smooth vari- 
ation in between, consistently with what was shown in lCappellari et alj (l2012h . The observed 
distribution of galaxy properties on the VP provides a clean and novel view for a number of 
previously reported trends, which constitute special two-dimensional projections of the more 
general four-dimensional parameters trends on the VP. We interpret the distribution of galaxy 
properties on the VP as due to a combination of two main effects: (i) an increase of the bulge 
fraction, with a corresponding increase of the concentration and decrease in Rc, which seems 
connected to the quenching of star formation, and (ii) dry merging, increasing galaxy mass 
and i?o by moving galaxies along lines of roughly constant ae (or steeper), while leaving the 
population unchanged. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: structure - galaxies: kinematics and dynamics 



1 INTRODUCTION 



Much of our understanding of galaxy formation and evolution 
comes from the study of dynamical scaling relatio ns relating galaxy 
luminosity or mass, si z e and kinema t ic (e. g. Faber & JacksonI 
19761; iKormendvl 1 19771: lOressler et alj 1 19871 : iFaberetalJ Il987l 



Diorgovski & Davislll987ir or regular trends in the distribution 



of g alaxy properti es as a function of t heir s c aling parameters 
(e.g. ISendereta l 1992; Burstein et"al] 1 19971 ; iKauffmann et alj 
l2003br i Gallazzi et al.. 2006), an d from the st udy of t heir evolution 

\\l996; Kelso n et alj[l997l ; 



with redshift (e.g. 



Ivan Dokkum et al. 



van D okkum & Fra__ , _ , 

1998l ; ITreu et al.ll2005i ; lFranx et allboOSh . 



The volume-limited ATLAS sample of nearby early-type 
galaxies jCappellari et alj[2011al . hereafter Paper I) constitute an 
ideal benchmark for studying scaling relations and the distribu- 
tion of galaxy properties, given the availability of a high-quality 
multi-wavelength dataset. In particular in ?, hereafter Paper XIX 
we used the state-of-the-art integral-field kinematics to derive accu- 
rate masses and global dynamical parameter. We found that galaxies 
lie, with very good accuracy, on a thin Virial Plane (VP) describing 
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galaxies in the parameter space defined by mass, velocity dispersion 
and projected half-light radius (AfjAM, o"e, ii™'"'). The existence 
of this plane is due to the virial equilibrium condition 



A/jam oc R" 



(1) 



and for this reason by itself it contains no useful information on 
galaxy formation. All the useful constraints on galaxy formation 
models come from the inhomogeneous distribution of galaxies in 
non-edge-on views of the VP and from the distribution of galaxy 
properties along the VP. 

This paper is devoted to a study of the non edge-on projections 
of the VP to see what we learn from it on galaxy format i on. Th is is 
done in the sp i rit of the classic papers bv lBender et"^ ( Il992h and 
iBurstein et However the fact that we have accurate dy- 

namical masses implies that our VP is extremely thin and it follows 
the virial equations quite accurately. For this reason we can ignore 
edge-on views of the plane and focus on non edge-on projections 
only. The thinness of the VP implies that any inclined projection 
show essentially the same information, after a change of coordi- 
nates. We can use standard and easy-to-understand observables as 
our main coordinates, instead of trying to observe the plane at a 
precisely face-on view. 
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In this paper, in Section 2 we summarize the sample and data, 
in Section 3 we present our projections of the VP. We illustrate the 
distribution of a number of quantities on the (M, Ue) and (A/, Rc) 
projection. We show the variation of the (M/L)jam, as well as 
of population indicators of MjL. We show the variation of galaxy 
concentration, intrinsic shape, morphology and stellar rotation. The 
variations of the IMF are separately presented in Section 4, together 
with a review of previous results on the IMF variation, in Section 5 
we discuss the implications of our findings for galaxy formation 
and briefly summarize our paper in Section 6. 



2 SAMPLE AND DATA 
2.1 Selection 

The galaxies studied in this work are the 260 early-type galax- 
ies which constitute the volume-limited and nearly mass-selected 
ATLAS'^'^ sample (Paper I). The object were morphologi cally se- 
lected as early-type accord i ng to the stan dard criterion JHubbld 
1 1 93 d ; Ide VaucouleursI 1 1 9591 : 1 S andagd 1 1 96 ll ) of not showing spiral 
arms or a disk-scale dust lane (when seen edge-on). The early- types 
are extracted from a parent sample of 871 galaxies of all morpho- 
logical t ypes brighter than Mk = —21.5 mag, using 2MASS pho- 
tometry dSkrutskie et al.ll2006l) . inside a local (D < 42 Mpc) vol- 
ume of 1.16 X 10^ Mpc^ (see full details in Paper I). 

In Paper XIX we compared our sample to previous samples 
for which accurate dynamical m asses have been determined ei- 
ther vi a gravitational lensi ng (e.g.lBolton et alJ|2008bl:lAuger et al 



or dynamics (e.g. iMagorrian et al.lll998l : Cappellari et al 



20061 ; i Thomas et a i] |2009h and conclude that it provides a major 



step forward in sample size and accuracy. 



2.2 Stellar kinematics and imaging 

Various multi-wavelengths datasets are available for the sample 
galaxies (see a su mmary in Paper I ). In this work we make 
use of the SAURON ( iBacon et an(200lh integral-field stellar kine- 
matics with in about one h a lf-ligh t radius Tie, which was in- 
troduced in 'Emsellem et al (E004), for the su bset of 48 early- 
types in the SAURON survey (de Zeeuw et alj [2002). The kine- 
matics of all galaxies in the ATLAS'^^ sample was homoge- 
neously extracted as described in Paper I, using the pPXF software 
Cappellari & Emselleml liooi) and the full MILES stella r library 



1"' 



Sanchez-Blazguez et alj 



2006l : lFalc6n-Barroso et alfcOIlh as tem- 



plates. 

Effective velocity dispersions CTc were measured by co-adding 
all the spectra in the SAURON datacube contained within the half- 
light isophote, or within the largest observed aperture. The result- 
ing "effective spectrum" was still fitted with pPXF using the full 
MILES stellar library as templates, and assuming a Gaussian line- 
of-sight velocity distribution. Gas emission around the possible H/3 
and [0777] emission lines were systematically excluded from the 
fits and the CLEAN keyword of pPXF was used to reject possible 
remaining outliers from the spectra. 

The photometry us ed in this work co mes from the Sloan Dig- 
ital Sky Survey (S PSS, pYork et aljlioool) data release eight (DR8 
lAiharaetalJl201l ') and was supplemented by our own photometry 
taken at the 2.5-m Isaac Newton Telescope in the same set of fil- 
ters and with comparable signal to noise for the rest of the sample 
galaxies dScott et alj|20ll hereafter Paper XXI). 



2.3 Measuring galaxy global parameters: MjL, Rc and ctc 

The dynamical masses used in this paper were obtained with 
the dynamical models of the ATLAS'^'^ sample presented in 
ICappellari et (12012) and described in more detail in Paper XIX. 
In brief the modelling approach starts by approximating the ob- 
served SDSS and INT r-band surface brightness distribution of the 
ATLAS galaxies^ using the Multi-G aussian Expansion (MGE) 
parametrization ( Emsellem et al ."1994'), with the fitting method and 
software of Cap pellari (2002Q. Full details of the approach and 
examples of the resulting MGE fits are given in Paper XXI. The 
MGE models are used as input for the Jeans Anisotropic MGE 
(JAM) modelling method^ which calculates a prediction of the line- 
of-sight second velocity moments (iifos) for ^ g iven set of model 
parameters and fits this to the observed Vru^„ dCappellaril I2OO8) 
using a Bayesian approach jGelman et a"i]|2004b . The models in- 
clude a spherical dark halo, which is parametrised according to six 
different sets of assum ptions (See Paper XIX for full details). In 
ICappellari et ^ i20I2h we showed that the adopted assumptions 
on the halo have insignificant influence on the measured trend of 
AI/L. For this reason in this paper we only use the two simplest 
sets of models from that paper (using the same notation): 

(A) Self-consistent axisymmetric JAM model, in which the dark 
matter is assumed to be proportional to the stellar one. We call the 
total (luminous plus dark) AI/L of this model (Af / L)jam; 

(B) A xisymmetric JAM model with spherical iNavarro et alj 
jl996l) (NEW) dark matter halo. Here the M/L of the stars alone 
(A//I/)stars is measured directly. 

iLablanche et alj j20I2l hereafter Paper XII) studied the accu- 
racy in the Af/L obtained with the JAM method, using A^-body 
simulations that resemble real galaxies, and concluded that the 
AI/L of unbarred galaxies can be measured with negligible bias 
< 1.5%, while the AI/L of barred galaxies that resemble typical 
bars found in the ATLAS^^ sample can bias the determination by 
up to 15% in our tests, depending on the position angle of the bars. 

The AI/L of the st ellar population ( M/L) saiD in the SDSS 
r-band was presented in ICappellari et It was extracted 

using the penalized pixel-fitti ng (pPXF) method and software of 
ICappellari & Emsellenj ("2004') for full spectrum fitting. In this case 
the model templates consisted of a regular r ectangular grid of 26 
ages and 6 metallicities [AI/H], assuming a ISalpeteil(ll955h IM F 
for reference, from the MILES models of IVazdekis et al. I I2OIOI) . 
We used the pPXF ke yword REGUL to enforce linear regulariza- 
tion (P ress et al .I199I eq. 18.5. 10) in the recovered set of templates 
weights from the fit. The regularization parameter was chosen for 
every galaxy to obtain an increase Ax^ ~ 1^/2 x A^pix in the x^, 
with respect to a non-regularized fit. In this way the recovered so- 
lution constitutes the s moothest one consist ent with the observed 
spectrum (see fig. 20 of lOnodera et alj|2012l for an example). This 
regularized approach helps reducing the noise in the (Af/L)saip 
determination, however the results are quite insensitive to the pre- 
cise value of this parameter and remain essentially unchanged for 
a range of plausible values. More details on our spectral fitting of 
the SAURON data and the resulting star formation histories will be 
presented in McDermid et al. (in preparation). 

Effective radii are defined as the projected half-light radii of 
the isophote containing half of the total analytic light of the MGE 
models from Paper XXI. Both circularized radii Rc and the more 
robust major axes J?™'*" were extracted as described in Paper XIX. 



Available from |http://purl.org/cappellari/idl] 
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M— a M-size 




0.35 0.48 0.61 0.74 0.87 



Figure 1. The Virial Plane and it projections. Tlie top two panels show the two main projections of the VP in the (Mjami ^^e) ™d (A^jaMi R^'^^) coordinates. 
Overlaid are hnes of constant (Te = 50, 100, 200, 300, 400, 500 km s~i (dashed blue), constant i?™^" = 0.1, 1, 10, 100 kpc (dot-dashed red) and constant 
= 10® , 10^ , 10^" , lO'^^ Mq kpc~^ (dotted black) predicted by the virial relation. The observed (Mjam , , R™'"') points follow the relation so closely 
that the coordinates provide a unique mapping on these diagram and one can reliably infer all characteristics of the galaxies from any individual projection. In 
each panel the galaxies are coloured according to the (LOESS smoothed) log(M/L)jAM, as shown in the colour bar at the bottom. Moreover in all panels the 
thick red Hne shows the ZOE relation given by equation (5), again projected according to the virial relation A/jam = 5.0 X (Tgi?™^''/G. While the top two 
panels contain different observable quantities, the bottom two panels merely apply a coordinate transformation to the quantities in the top two panels, to show 
the effective phase space density /off = l/(crRc'^) and effective mass surface density Se = Af jAM/(27ri?c^). Two galaxies stand out for being significantly 
above the ZOE in the (Mjam, fe) and (AfjAM, Se) projections. The top one is NGC 5845 and the bottom one is NGC 4342. 



3 PROJECTIONS OF THE VIRIAL PLANE 
3.1 Total M /L variations 

We have shown in Paper XIX that the existence of the FP is almost 
entirely due, with good accuracy, to a virial equilibrium condition 
combined with a smooth variation in M/L. Once this is clarified, 
the edge-on projection of the Virial Plane becomes uninteresting 
from the point of view of the study of galaxy formation, as it merely 
states an equilibrium condition satisfied by galaxies and it does not 
encode any memory of the formation process its elf. This is in agree- 
ment with previous finding s with simulations ( lNiDotietal.ll2003l : 
IBovI an-Kolchin et al.ll200q) . All information provided by scaling 
relations on galaxy formation is now encoded in the non edge-on 
projections of the Virial Plane, and first of all in the distribution 
of M/L on that plane. In Paper XIX we also confirmed that M/L 
correlates remarkably tightly with cje (Cappellari et al. 2006). This 
is especially true (i) for slow rotators, (ii) for galaxies in clusters 
and (iii) at the high-end of the Oe range. Here we look at the entire 



Virial Plane and try to clarify the reason for these and other galaxy 
correlations. 

In a classic paper lBender et all ( Il992h studied the distribution 
of hot stellar systems in a tiiree-dimensional space, they called k 
space, defined in such a way that one of the axes was empirically de- 
fined to lie nearly orthogonal to the plane. This made it easy to look 
at both the edge-on and face-on versions of the plane. In this paper, 
thanks to the availability of state-of-the-art integral-field kinemat- 
ics and the construction of detailed dynamical models, we can use 
mass as one of the three variables (A/jam, o"e, Re)- We have shown 
that in these variables the plane is extremely thin and follows the 
scalar virial equation Mjam = 5.0 x J?™^''/G within our tight 
errors. This implies that any projection of the plane contains the 
same amount of information, except for a change of coordinates. 
Instead of looking at the plane precisely face-on, we decided to con- 
struct special projections that correspond to physically-meaningful 
and easy-to-interpret quantities. 
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Our selection of meaningful projections of the VP is shown in 
Fig.[T] We use as horizontal axis in all plots our main mass variable 



■c,b 



MjAM = LX (M/L) JAM ~ 2 X Ml/2 ~ 



(2) 



where (M/L) jam is the total (luminous plus dark) dynamical 
M/L obtained using self-consistent JAM models, L is the total 
galaxy luminosity and M1/2 is the total mass within a sphere of 
radius rij2 enclosing half of the total gala xy light, w here r-1/2 ~ 
1.33i?o l iHem quist 1990; Ci otti 1991; W olf et alf201(]l Paper XIX). 



Although the self-consistency assumption, where the total mass is 
proportional to the stellar mass, is not justified at large radii, it is 
accurate within the re gion where we have st ellar kinematics (about 
IRc)- It was shown in lWilliams et and in Paper XIX that 

(M/L) JAM closely reproduces the total (M/L)(r = Re) inside 
a sphere (actually an iso-surface) with mean radius the projected 
half-light radius Rc, derived using models which explicitly include 
dark matter. Different assumptions on the dark halo produce minor 
differences in {M/L){r = R^). Given that {M/L){r = Ra) is 
nearly insensitive to the choice of the halo assumptions, we choose 

the self-consistent one, being the simplest. 

As illustrated in Paper XIX (see also IWilliams et ai]|2009l) , 
most of the galaxies in our sample are consistent with having 
small fractions of dark matter within a sphere of radius r — Rc, 
with a median value of just 12% of DM within that radius for 
our standard models (B). This implies that M1/2 is dominated by 
the stellar mass. For this reason Mjam is a quantity that very 
closely represents and is directly comparable to the total stellar 
galaxy mass used in numerous previous studies. Stellar mass seems 
to relate well with g alaxy properties and is of ten used to study 
galaxy formation (e.g. lKauffmann et aI.ll2003al lbl; lHvde & Bernardil 
l2009bl) . The difference of our mass parameter is that it does 
not suffer from th e uncertainties related to the s tellar popu- 
lation models (e.g . iMaraston et aO |2006|; Gallazzi & Bell 200' 
IConrov et alj|2009l ; iLonghetti & Saraccoll2"o09i ; iWuvts et alj(200' 



more over it automatically includes the effects of a non-u niversal 
IMF jvan Dokkum & Coru-ovll2010l ; ICapDeIlari et al.ll2012l) . Being 
a measure of the total enclosed mass within a spherical region, and 
thus being directly related the dynamics, Mjam represents the ideal 
parameter which one would like to use in scaling relations. Note 
that AfjAM, unlike mass determinations obtained via strong lens- 
ing, does not include the possible contr ibution of dark matte r along 
the cylinder parallel to the line-of-sight jDutton et al.l2011al) , which 
provides an useful additional constraint on dark matter at large radii, 
but complicates the interpretation of scaling relations in the galaxy 
centres. 

In the top-left and top-right panels of Fig. [T] we show the pro- 
jections of the VP along the (Mjam, o"e) axes. The colour in these 
diagrams represent the dynamical M/L inside a sphere of radius 
_Rc. Three important results are clear from these plots: 

(i) Both projections are equivalent and provide basically the same 
picture, apart from a coordinate transformation. This is expected 
from the tightness of the VP. Moreover, given that the VP nearly 
follows the virial relation, one can plot virial coordinates defined 
by Mjam = 5.0 x a^Rc/G, to construct a consistent and accu- 
rate system of coordinates in both panels to infer other dynamical 
quantities; 

(ii) Galaxies define a clear zone- of-e xclusion (ZOE), as al ready dis- 
covered bv lBender et alj ( Il992h and lBurstein et al. ( Il997l) . however 



we find a clear break at a mass M ~ 3 x 10 Mq and two nearly 
power-law regimes above or below this value. The ZOE is approxi- 
mated by the equation 



/ Mjam ' 




'1 


-I 


( Mjam \ 


\ MjAM.b ^ 




2^ 




\MjAM,\> ) 



-, (/3-7)/q 



(3) 



0.75 



7 

MjVm above MjAM.b 



-0.30. The re- 



cx m; 



below 



with i?c,b = 0.7 kpc, a = 8, /3 
lation has an asymptotic trend R^ (x 
2.8 X 10^'^ A/0, and a sharp transition into R^ 
this break. The values were determined by simultaneously match- 
ing the observed boundary in the galaxy distribution in all VP pro- 
jections. These values are close to those already reported in Paper I 
(eq. 4 there), using Jf-band luminosity in place of mass and 2MASS 
i?c instead of MGE ones. The maximum in the galaxy density at 
MjAM,b that we infer from our sample is also clearly vi sible in 
the much larger SDSS sample of Ivan Dokkum et ani2008l fig. 2). 
This shows that the relation is robust. The ZOE relation can be con- 
verted into a (MjAM,o"e) one, or into other projections using the 
scalar virial relation AfjAM = 5.0 x a^R^^^ /G. The location of 
the break we find agrees with the value at which scaling relations 
of large sample of SDSS galaxies show a subtle deviation from a 
straight line (Hyde & Bernardi 2009a). Interestingly in both cases 
th e characteristic mass at the cusp coincides with the value reported 
bv lKauffmann et alj ( l2003bl) . as the fundamental dividing line be- 
tween the two distinct families of passive and star-forming galaxies, 
(iii) The contours of constant (M/L) jam closely follow lines of 
constant above ^ 100 km s~^. A comparable agreement is 
reached using lines parallel to the ZOE i?f oc A/jVm. Neither 
mass A/jam, nor size Rc or surface mass density Ee provide a 
comparably good approximation to the {M/L)jam contours, al- 
though Se provides a better approxima tion than the other two. This 
is consistent and explains the finding bv lCappellari et all j2006l) that 
(Je and not mass or luminosity, is the best tracer of (M/L) jam. 
At lower ae < 100 km s"^ the (M/L) jam contours start deviating 
from the lines of constant a and tend to lie closer to lines of con- 
stant i\/jAM. Given that the A//L describes the deviations between 
the FP and the VP, this twist in the contours demonstrates that the 
FP is warped with respect to the VP, and explains the sensitivity of 
the FP parameters to the region included in the fit Ip'Onofrio et alj 
.2008 ; Gargiulo et al. 2009; Hyde & Bernardi 2009^^ 

Two galaxies stand out for being significantly above the ZOE 
in the (MjAM,o"e) and (AfjAM,Se) projections. The top one is 
NGC 5845 and the bottom one is NGC 4342. The two objects have 
very high surface brightness and consequently excellent kinematic 
and photometric data. They are genuine examples of dense ob- 
jects in t he nearby Universe . These two outliers were already pre- 
sented in ICaDpeIlaril|2011ah and their compactness was later dis- 
cussed bv IJiang et alj ((2012). However they have smaller masses 
(A/jam ~ 3 X 10^° A/p)) than their high -redshift counterparts 
jCimatti et alj2008l ; lvan Dokkum et al.l2008l) . In all panels of Fig.[I] 
the measured (M/ L) jam have been adaptively smoothed using the 
Loceilly Weighted Regression robust technique (dubbed LOESS) of 
ICIevela na ('1979'), straight forwardly generalized to two dimension 



Cleveland & DevliiJ ( fl98a) . In all plots we adopt a regular- 



ization factor / = 0.3, and a linear local approximation, which 
implies that a fraction 30% of the points are adopted every time 
for the tricube-weighted robust fit of a local plane. To deal with 
the different scales of the axes, and the elongated relations, before 
applying LOESS we re-normalize the coordinates so that their ellip- 
soid of inertia reduces to a circle. We do not show here the original 
data, as the scatter in (M/L)jam is already analysed in detail in 
Paper XIX. Given that most of the variation in the smooth surface 
happens orthogonally to the constant ae lines, the scatter from the 
smooth surface of 27% is nearly unchanged from the value of 29% 
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in the global [M/ L) — tJe relation, further confirming that all the 
variation happens orthogonally to cTe- 

In the bottom left and right panels of F ig.[T]we show the effec - 
tive phase-space density, defined following iHemquist et alj ( Il993h 
as /eff = l/((Te-Rc^) and the mass surface density, defined as 
Ee = MjAM/(27ri?e^). Notc that, while we use i?™''" for the 
mass-size plane, one has to use Rc for Eg. The two bottom pan- 
els do not plot new data, but are obtained by rearranging the three 
variables shown in the top panels, however they illustrate how 
the above trends relate to other physical quantities. In particular 
the phase-space density is interesting because it can only decrease 
during collisionless galaxy mergers, due to Liouville's theorem 
l lCarlberg|[T98^ : iHernquist et Zlll993h . Interestingly the ZOE we 
find is nearly flat in /eff below MjAM.b = 3 x 10^" (a value 
7=1/3 gives /off = const in the virial case) and starts decreas- 
ing at larger masses, as one would expect when dry mergers start 
becoming more important. 



3.2 Stellar population indicators of {M / I/)pop 

We have shown in Paper XIX that the stellar matter dominates in 
the regions we study. If this is indeed the case we would expect 
stellar population indicators of the (i\f/L)pop of the stellar popula- 
tion to closely follow the behaviour of the dynamical (M/L)jam. 
Before addressing this question with our own data, the reader is 
st rongly encourag ed to compare Fig.[T]to the right panel of figure 15 
of lGallazzi et alj (2006), which shows the luminosity-weighted age 
versus the stellar mass in the (a, hh) projection. Even though 
Fig. [T] uses dynami cal quantities, which are measured via dynam- 
ical models, while iGallazzi et alj ( l2006h derives population from 
line strengths, the two plots resemble each other very closely, in- 
cluding the change of orientation of the contours around a ~ 100 
kms~^. This comparison already indicates a close link between 
M/L and stellar popul ation (mainly age), in broad agreement with 
ICappellari et"aLl l l2006h . Impo rtantly, given that our sample is 100 
times smaller than the one of lGallazzi etai] (|200^, this compari- 
son also confirms the ability of the LOESS technique to recover the 
underlying distribution from our smaller galaxy sample. 

We now address the relation between dynamical and popu- 
lation M/L using our own dataset and a different approach. Two 
main simple tracers of (M/L)pop have been proposed in the past: 
(i) the B — R colour, which was shown by (Bell & de Jong 2001) 
to trace (M/L)pop alone, for a wide range of metallicities, and (ii) 
the line-strength index Hj3, which was shown to satisfy a simi- 
lar property in Cappellari et al. (2006) for a wide range of metal- 
licities (see also iWorthevI 1 1 994 ) . Here we adopt the HjS determi- 
nation for our ATLAS^'' galaxies, as derived from the very high 
S/N ef fective SAURON spect ra. The parameters for the SAURON 
survey (Ae, Zeeuw et al. "2002*) subset of 48 galaxies was already 
given in Kuntschner et al. ( 2006). The homogeneous extraction for 
the ATLAS^'^ sample will be given in McDermid et al. (in prepa- 
ration). As choice of colour we use the SDSS g — i one, which is 
av ailable for 223 ATLAS^ ° galaxies in SDSS DR8 and was shown 
bv lGallazzi& Belli ( l2009h to provide on average the smallest uncer- 
tainties and the most stable results, among the SDSS bands. 

The results are presented in Fig. |2| and show a quite good 
agreement between simple population estimators of M/L and the 
dynamical one. Both HP and colour are nearly constant along lines 
of constant a as it was the case for Fig. [Tj This confirms the fact 
that the key driver of the total M/L is the stellar population. We 
do notice some systematic differences at the low-cr end, with the 
contours of constant (i\f/L)pop not quite following the popula- 



tion trends. As the effect is visible both in colour and Hji we be- 
lieve it is significant. It may be due to the presence of extended 
star formatio n episodes, which a r e expected to be m ore prevalent 
at low mass faeavens et ai]| 20041 : Ifhomas et alj|2005l) and cause a 
break in the relation between population observables and M/L. The 
study of the relation between our parameters and the actual phys- 
ical parameters of the population, like age, metallicity and abun- 
dance ratios, is presented in McDermid et al. (in preparation). Anal- 
ysis of the stellar pop ulation parameters for the SAURON survey 
fde Zeeuw et al. 2002) subset of 48 galaxies was already presented 
in C Kuntschner et aU201Qi) . 



3.3 Variations of galaxy concentration and morphology 

In the previous sections we showed that galaxy properties follow 
smooth trends in the projections of the VP. We have seen that the 
dynamical M/L near the galaxy centre (IRc) is driven by a varia- 
tion of the stellar population, with the dynamical M/ L constituting 
an accurate tracer of the galaxy population! These trends distribute 
galaxies in the VP along lines of nearly constant de. In Paper I 
we showed that the luminosity-size relation changes gradually with 
sizes decreasing as a function of the bulge ratio, indicated by the 
morphological classification. We interpreted the apparent size de- 
crease as an actual increase of the bulge mass between different 
morphological types. In Paper II and III, we showed that the vast 
majority of local ETGs is constituted by disk-like systems, the fast 
rotators, and in jCappellari et alj|201 lb . hereafter Paper VII) we 
discussed how fast-rotators morphologically resemble spirals with 
a variety of bulge fractions. We concluded that fast rotators simply 
constitute the end point of a smooth sequence of disks with increas- 
ing bulges, and this lead us to propose a fundamental change in 
our view of galaxy morphology (Paper VII). If this picture is cor- 
rect, we should be able to find an indication for a systematic change 
in the bulge fraction, in the fast rotators class, while moving from 
the region of the VP populated by spiral galaxies towards the ZOE, 
dominated by the oldest and reddest galaxies, with the largest M/L. 

A detailed photometric bulge-disk decomposition of the un- 
barred galaxies subs et of the ATLAS'"^ sample is presented 
iKrainovic et alj j2012l hereafter Paper XVII). Although it provides 
a well established definition of the bulge luminosity, it is a complex 
technique which depends on the extraction details and may suffer 
from degeneracies. Moreover it is unclear how it should be applied 
to barred galaxies. For this reason here we look instea d at the sim- 
ple concentration parameter TGC dTruiillo et al.ll200l]) . It measures 
the ratio between the flux within 1 J?e and 7?e/3 which we extracted 
from the MGE photometric models as described in Paper XIX. For 
a given galaxy mass an increase of the bulge fraction produces by 
definition a larger concentration. The distribution of galaxy concen- 
trations on the VP is shown in the top panel of Fig. [3] Even though 
the values are noisier than the clean M/L trends, a systematic vari- 
ation is evident, with galaxy concentration roughly following the 
M/L, Hp and g ~ i colour trends, which all tend to be constant 
along lines of constant Ue. This trend is cons istent with the trends o f 
concentration against Ee and tJe reported bv^ Graham et alj ( f2001bl) . 
and the corres ponding correlation with supermassive black holes 
^Graham et aLil2001iJ). It is also consistent with a similar trend be- 
tween M/L and lSersi^ j 1 9681) index reported bv lD'Onofrio et al] 
j20I ih . Here we interpret all these trends as due to the increase of 
the bulge fraction, rather than a smooth change in the profile shape 
of a spheroidal system. It's important to stress that our term "bulge" 
refers to an increase of the mass fraction in the centre of the galaxy. 
It does not necessarily imply a correspondence with the definition 
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Figure 2. Indicators of (M/L)pop on tiie Virial Plane. Top Panels: Same as in Fig.[T] with colours showing the log Hfi line-strength, from McDermid et al. 
(in preparation), which is a good indicator of (M/L)p op. The open circles represent a de termination for galax ies at redshift z > 1.4: blue symbols from 
Cappellari et al.l i2009l) . green symbol fromlCenarro & T ruiillo (2009), magenta symbol from lOnodera et al]j2012l) and black symbol from jvan de Sande et all 
201 Not shown is the cr value from lvan Dokkum et alj (2009), which is off the vertical scale. Bottom Panels: Same as in Fig.[T] with colours indicating the 
SDSS galaxy colour g — i. Note that both (M/L)pop and g — i show the same trends as (Af/L)jAM in Fig.[T] 



of "bulge" that one infers via classic bulge-disk decomposition (Pa- 
per XVII), namely a light excess over an outer exponential com- 
ponent. There is some evidence that the concentration may follow 
population even more closely than M / L, possibly due to the con- 
tribution of extended star formation episodes as we commented for 
Fig. [2] 



The concentration plot demonstrates that indeed, even within 
the ETG class (excluding the few slow rotators), the variation 
of galaxy properties follows the bulge growth. Interestingly the 
plot also shows a slight decrease in the concentration at the 
largest mass end Mjam^2 x 10^^ M©. This value coincides 
with another charact eristic scale in ETGs, above which galaxies 
start to be rounder jTremblav & Merrlt^ 1 19961 : Ivan der Wei et alj 
l2009bl). to have flat central profiles in their su rface brightness 
jpaber et ah , 1997: Graham 2004: Fe rrarese et al.i i2 006) and to de- 
viate from power-law colo ur-magnitude relations jperrarese et alj 
l2006l : lBernardi et al.ll201 ll) . 



3.4 Variations of galaxy intrinsic flattening and rotation 

Our dynamical modelling effort provides the unique opportunity 
of having the galaxy inclination for the entire ATLAS^'^ sample 
of galaxies. We verified that the JAM inclination agrees with the 
one inferred from the geom etry of the dust, for 26 galaxies with 
regular dusts disks (see also ICapDellariii2008.) . We further verified 
that we can recover with JAM models the inclination of simulated 
galaxies that resemble typical objects in our sample (Paper XII). 
Although one should not expect the inclination to be reliable for 
every individual galaxy, it is expected to be accurate for most of 
our sample. The inclination allows us to recover the intrinsic shape 
of individual objects. The deprojection is done using the elliptic- 
ity of the outer isophote, measured in Paper II at radii typically 
around 4_Rc, at the depth limit of the SDSS photometry. We then 
computed the intrinsic ellipticity ass uming an oblate spheroid ge- 
ometry as teinnev & Tremairi3l2008l) 



£intr = 1 - Vl + e(e-2)/sin2i, (4) 
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Figure 3. Concentration and shape on the Virial Plane. Top Panels: Same as in Fig.[T] with colours indicating the concentration parameter TGC of lTruiillo et alj 
jloOl). The trends are noisier than M/ L but there is a clear increase in the concentration along lines of roughly constant ae and a slight decrease above a 
characteristic mass M}A_M > 2 X 10" Mq . Middle Panels: Same as in Fig.0 with colours indicating the intrinsic ellipticity Cintr of the galaxies outer regions, 
computed by de-projecting the observed ellipticity e at large radii (Paper II) using the inclination inferred from the JAM models. Except again for the special 
region above Mjam ^ 2 x 10^^ Mq , all the rest of the ETGs have on average the flattening of disks e 0.65 consistently with the Monte Carlo approach of 
Paper III, which suggests fast rotators have a typical intrinsic ellipticity eintr ~ 0.7 and the statistical inversion in Weijmans et al. (in preparation). 



© 2012 RAS, MNRAS 000.[T]l25] 



The ATLAS^^ project - XX. Mass-size 9 



M—a M— size 




0.13 0.26 0.38 0.51 0.64 

M— CT M— size 




0.05 0.07 0.10 0.12 0.15 



Figure 4. Same as in Fig.[T] witli colours indicating tiie specific angular momentum parameter A/j of iEmsellem et al.l2007l) . as given in Paper III. Top Panels: 
all ATLAS^°galaxies are included. Bottom Panels: Only the slow rotators are shown. They seems to define a sequence in the VP, and they show a decrease in 
the rotation above A/jam ^ 2 x IO^^Mq. 



The results of the deprojection is shown in the middle panels of 
Fig- El Contrary to all previous diagrams the distribution of the in- 
trinsic ellipticity for ETGs on the VP shows a completely differ- 
ent trend. Few galaxies stand out in the top-right comer for being 
nearly round, with intrinsic ellipticity Sintr ~ 0.2. These galax- 
ies are all located above the same characteristic mass Mjam ^ 2 x 
W^^A4q, where the top panels showed a drop in the central den- 
sity. This result confirms previous statistical studies of ETGs shapes 
jTremblav & Merritlll99a : lvan der Wei et alj2009bl : lBemardi et alj 
I20T l|). Below this characteristic mass, the mean ellipticity drops 
dramatically and sharply to about eintr ~ 0.65, characteristics of 
disks. This galaxy-by-galaxy shape deprojection confirms a similar 
result on their shape distribution inferred via Monte Carlo simu- 
lations from the distribution of fast rotators ETGs on the {Xr, e) 
diagram, suggesting fast rotators have a typical intrinsic ellipticity 
Eiiitr ~ 0.7 (fig. 15 in Paper III). It also agrees with a statistical in- 
version of their shape (Weijmans et al. in preparation). This agree- 
ment also seems to confirm the reliability of the inclinations derived 
via the JAM models. The result strongly confirms our previous con- 
clusions that fast rotators as a class are disk like (Paper II, III, VII). 
The fact that the outer parts of fast rotators are as flat as disks, while 



the central ones show an increase of concentration, conclusively 
demonstrates that the latter is due to an increase of the bulge (or 
spheroid) mass, and not to the smooth variation of the profile in a 
spheroidal object which follows a R^^"^ profile. At the lowest mass 
range, near the ZOE, there seems to be a marginal decrease of eintr- 
The significance or reality of this feature is however unclear. 

The round and massive objects that stand out from the dis- 
tribution are the prototypical slow rotators, which also stand out 
for their slow specific stellar angular momentum, as shown in Pa- 
per III. Another view of that fact is shown in the bottom panel 
of Fig. [3] which plot s the angular momentum parameter Xr of 
lEmsellem et alj | |2007|) for the slow rotator ETGs only, as tabulated 
for the galaxies of the ATLAS^'^ sample in Paper III. The slow ro- 
tators are found to populate a rather narrow sequence on the VP, 
with a small range of Ee. There is a clear trend of decreasing rota- 
tion with increasing mass as noted in Paper III, or possibly a sharp 
transition around Mjam ~ IO^^Mq, which is similar to the char- 
acteristic mass defined by concentration and intrinsic ellipticity. Al- 
though we tend to cover a smaller field with our kinematics for the 
largest galaxies, the qualitative difference in the kinematics for the 
most massive slow rotators is already clearly visible from the ve- 
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locity fields (Paper II) and it is not due to tiie difference in the field 
coverage. 

In summary the distribution of galaxy properties on the 
(MjAMjfe) and (MjAM,i?™^^) projections of the Virial Plane 
can be understood as due to a smooth variation in the bulge mass 
as already suggested in figure 4 of Paper I. The connection between 
bulge mass and galaxy properties, and the close link between ETGs 
and spiral galaxies is further illustrated in Fig. [S] which includes 
the location of the spiral galaxies of the ATLAS^° parent sample 
(Paper I) together with the ETGs. The masses of spiral galaxies was 
approximately estimated assuming a fixed M/Lk ~ 0.8 MqILq, 
which ensures agreement between Mjam and the JC-band lumi- 
nosity at the lowest masses. The plot shows that ETGs properties 
vary smoothly with those of spiral galaxies. Galaxies with negligi- 
ble bulges are almost invariably star forming and classified as late 
spirals. Galaxies with intermediate bulges can still form stars and 
be classified as early-spirals, or can be fast rotator ETGs. But the 
galaxies with the most massive bulges, as indicated by their largest 
concentration or (Je are invariably ETGs. They have the largest 
M / L (Fig.[T] the reddest colours and smallest H/? (Fig. |2l but are 
still flat in their outer parts, indicating they have disks (Fig. [3) and 
generally still rotate fast (Fig. |4j. An exception are galaxies with 
masses Mjam ^ 2 x lO^^M©, which stand out for being nearly 
round, and non rotating. 



4 SYSTEMATIC VARIATION OF THE IMF 

Although the IMF variation could have been included among the 
other galaxy observables described in Section|3] we keep the IMF in 
a separate section for a more detailed coverage than for the other ob- 
servables. In recent months there has been a large amount of interest 
on the IMF variation. Before describing our new ATLAS^° results, 
in the following section we provide a summary of this rapidly evolv- 
ing field. To clarify the significance and robustness of the various 
recent claims, including our own, we place particular emphasis on 
the modelling assumptions that went into the various studies. 

4.1 Summary of previous IMF results 

4.1.1 IMF of the Milky Way 

The stellar initial mass function (IMF) describes the distribution of 
stellar masses when the population formed. Nearly every aspect of 
galaxy formation studies requires a choice of the I MF to calculate 
predictions for galaxy observables ( lKennicuttll998h . For this reason 
the IMF has been subjected to numerous investigations since the 
first d etermination m ore than half a century ago finding it has the 
form ( ISalpeteilll955h 

("(m) oc = m~^'^ , (5) 

for m > 0.4Mq, where m is the stellar mass. 

In the Milky Way the IMF can be measured via direct stel- 
lar counts. Various determination in different environments lead to 
the finding of a remarkable universality in the shape of the IMF 
( lKrouDall2002l) , with the IMF being well described by the Salpeter 
power slope x = —2.3 for m < CSAf© and a shallower on e 
at smaller masses x = — 1.3 for m > O.SMq ( iKroupallioOlh . 
The Milky Way IMF can also be described by a log-normal dis- 
tribution ( Chabrier 2003), which has some theoretical justification. 
However the Kroupa and Chabrier IMFs are essentially indistin - 
guishable from an observational point of view jKroupa et alj20l3) . 
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Figure 5. The Mass-size distribution for spiral galaxies and ETGs. The 
ETGs of the ATLAS^^ sample are coloured according to their (M/L)jam 
as in Fig. [T] This plot however also includes the approximate location 
of spiral galaxies in the ATLAS^'^ parent sample, where masses are ob- 
tained from i^T-band luminosities, scaled to match A/jam (2 plot with both 
ETGs and spirals luminosities was shown in Paper I). Early spirals (Sa-Sb: 
T < 4) are indicated by the small black filled circles, while late spirals 
(Sc-Irr: T > 4) are shown with the small filled red circles. Late spirals are 
larger than ETGs, while early spirals overlap with the ETGs with low M/L. 
There are essentially no spirals in the region of the oldest and reddest ETGs, 
which have the largest M/L. This illustrates the fact that the increase of the 
spheroid is required to make a galaxy old and red and consequently produce 
the largest M /L. 

Recent IMF results were reviewed bv lBastian et alj ( l2010h who con- 
cluded no clear evidence existed for a non-universal IMF in our 
galaxy and among different galaxies. 

4.1.2 IMF from ionized ga.s emission or redshift evolution 

In external galaxies individual stars cannot yet be resolved down to 
sufficiently small masses for IMF studies. Some indirect constraints 
on the slope of the IMF above ~ IA/0 can be obtained by combin- 
ing observations of the Hq equivalent widths, which is related to 
the number of ionizing photons, with the galaxy colour, w hich is 
a function of t he galaxy stellar populatio n jKennicutlll 19981) . Using 
this technique iHoversten & Glazebrookl ( 1200 8 l) inferred a variation 
of the IMF with more massive galaxies having a more top-heavy 
IMF than Salpeter, a result confirmed under the same assumptions 
jGunawardhana et alj 1201 iT) by the GAMA survey jPriver et alj 
l201lh . 

Constraints on the IMF can also be obtained by comparing the 
local (z = 0) stellar mass density to the integral of the cosmic 
star formation hi story. Iwilkins et"al] ( l2008h finds that the local stel- 
lar mass density is lower than the value obtained from integrating 
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the cosmic star formation history (SFH), assuming that all the stars 
formed with a Salpeter IMF. They propose a ti me- variable IMF to 
reconcile the observations (see also lDavj|2008l) . 

Alternative constraints on the IMF can be obtained by study- 
ing the evolution of the M/L normalization in samples of galax- 
ies as a function of redshift, for example from the evolution of 
the tilt and normalizati on of the Fundamental Plane with redshift 
( lRenzini&Ciottill993i) . From these constraints and from the metal 
abundance in the clusters interstellar medium [Renzinil ( |2005|) con- 
cluded that elliptical galaxies as a class must have an IMF close to 
Salpeter for stellar masses m > IM©, while the IMF must flatten 
at lower masses, similar to the Kroupa IMF 

A related approach was employed bv lvanDokkuml ( l2008 l) who 
used the ratio of luminosity evolution to colour evolution of massive 
galaxies in clusters to constrain the IMF. He originally concluded 
that the IMF must be top heavy at high redshift, but this result was 
recently revised to conclude that the IMF a bove m > IMp) is actu- 
ally n ot inconsistent with a Salpeter slope dvan Dokkum & Conrovl 
|2012|) . 



4. 1.3 tensing or dynamics IMF constraints for external spirals 

A more direct approach to constraining the variation of the IMF 
in galaxies consists of measuring the mass of the stars, and com- 
paring this with the predictions of stellar population models. It is 
important to stress however that this method does not measure the 
shape of the IMF directly, but only its overall mass normalization. 
The method can only verify whether the stellar mass in a certain 
region within a galaxy is consistent with that due to a stellar pop- 
ulation with a certain assumed IMF. Moreover these methods do 
not measure the same IMF one can infer via direct stellar counts in 
star clusters. They measure instead the stellar mass distribution due 
to the superposition of the IMFs from a large number of star clus- 
ters in the galaxy as well as the IMFs coming from the accretion of 
satellite galaxies during the galaxy hierarchical growth. As pointed 
out by Kroupa & Weidner (2003), even if the IMF was universal 
in every cluster in a galaxy, the integrated galactic IMF (IGIMF) 
would be different from the universa l one, due to the no n-uniform 
mass distribution of star clusters (see lKroupa et al.l|20I2l . for an in- 
depth discussion). Hierarchical evolution complicates the picture 
even further For these reason, when we state e. g. that "a galaxy 
is consistent with the Salpeter IMF", we simply imply, consistently 
with all previous similar studies, that "the galaxy has the same stel- 
lar M/L of a stellar population with the given age and metallicity 
(and abundances), and the Salpeter IMF". 

The first convincing constraint on the IMF of external galaxies 
was obtained for a sample of 21 spira l galaxies, using the kine- 
matics of the gas. iBell & de Jon3 ( 1200 1[) concluded that if the IMF 
is universal, a sensible assumption at the time, it cannot have the 
Salpeter form, but it must be lighter, consistently with Kroupa-type. 
The result was later confirmed, sti ll using gas kinemat ics, for a sam- 
ple of 34 bright spiral galaxies bv lKassin et al.l ( |2006|) . 

Th need for a light IMF was inferred for the Einstein Cross 
spiral galaxy gravitation al le ns system (H uchra et al. 1985) by both 
Ivan de Ven et al. andlFerreras et alj ( i201Q) . and for another 

spiral le ns galaxy by Suvu et alj ( l2012h ~As part of the DiskMass 
Survey dSershadv et alj 20IOl) . which obtained integral-field stel- 
lar kinematics of a sample of 30 spiral galaxies, a light IMF 
was indirectly confir med from the sub-maximality of the disks 
jSershady et alj2()Tl ). As part of the SWELLS survey of spiral lens 
galaxies i Treu et alj|201 ll) the need for a light IMF for most spiral 



gala xies was also found from the analysis of 20 strong gravitational 
lens dBrewer et aTll2012h . 

An inconsistency of the Salpete r IMF norma l ization for low 
mass galaxies was also inferred by iDutton et alj ( 1201 Ibl) , using 
simple galaxy models with a bulge and a disk, trying to repro- 
duce global trends for a large sample of galaxies extracted from 
the SDSS. 

In summary there is a good agreement on the fact that spiral 
galaxies as a class must have a normalization lighter than Salpeter 
and similar to Kroupa/Chabrier. This is a robust result due to the 
fact that a Salpeter IMF would over-predict the total mass in the 
galaxy centres. It is unclear whether some spiral galaxies have a 
Salpeter IMFs, and of course wh ether the IMF varies within the 
galaxies themselves, as claimed bv lDutton et alj ( 1201 2cl) . 

4. 1.4 Lensing or dynamics IMF constraints on early-types 

An early attempt at constraining the IMF of 21 elliptical galaxies, 
using detailed spherical dynamical models including dark matter, 
found a general consistency between the stellar (i\f/L)stars and the 
one of the population (M/L)pop, using the Kroupa IMF, but could 
not accurately constrain the IMF n ormalization due to large obser- 
vational errors ( lGerhardetal .' 2001). A large study of SDSS ellip- 
tical galaxi es using fixed spherical H ernquist galaxy models with 
dark halos dPadmanabhan et al]|2004l) found a mass excess over the 
predictions of stellar population models with a fixed IMF, increas- 
ing with luminosity. This was interpreted as a incre ase of the dark 
matter fraction. A similar conclusion was reached bv lZaritskv et alj 
( l2006h while studying the fundamental manifold. A caveat of those 
studies was the use of homologous stellar profiles or approximate 
assumptions to study systematic variations in the heterogeneous 
galaxy population. 

To overcome thes e limitations, as part of the SAURON project 
dde Zeeuw et alj[2002l) . we constructed self-consistent axisymmet- 
ric models reproducing in detail both the photomet ry and the state- 
of-the -art integral-field stellar kinematics data ( Ems ellem et alj 
|2004|) for 25 early-type galaxies (Cappellari et al. 200^. The as- 
sumption that mass follows light is not accurate at large radii, but 
it is a good assumption for the region where kinematics is available 
( < IJie), and provides accurate measurements of the total (lumi- 
nous plus dark) M/L within a sphere of radius r ~ (see Pa- 
per XIX). The resulting improvement in the M/L accuracy and the 
removal of systematic biases allowed us to strongly confirm that 
"the total and stellar M/L clearly do not follow a one-to-one rela- 
tion. Dark matter is need ed to explain the diffe rences in M/L (if 
the IMF is not varying)" dCappellari etai] |2006'. see fig. 17 there). 
Although dark matter seemed, at the time, still a more natural expla- 
nation of our observations, if the measured trends are re-interpreted 
as an IMF variation, they would imply an IMF heavier than Salpeter 
for some of the oldest objects. That study also concluded that ;/ 
the IMF is universal it must have a mass normalization as low as 
Kroupa IMF, consistently with the results for spiral galaxies, other- 
wise the stellar mass would overpredict the t otal one for a number 
of galaxies. A similar finding was obtained bv lFerreras e7ZI ( l2008l) 
using gravitational lensin g of 9 elliptical galaxies extracted from 
the SLACS survey dBolton et al.l2006h . 

A number of subsequent studies used approximate models 
reproducing the velocity dispersion of large galax y samples to 
study dark matter and IMF variations in galaxies. Tortora et alj 
d2009l) used spherical isotropic models with Sersic (1968) pro- 
files to reproduce the large velocity dispersion compilation by 
IPrugniel & SimienI ( fl99^ . I Graves & Fabej d2010h used the Funda- 
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menta l Plane of early-type galaxies from SDSS data. ISchulz et alj 
combined weak lensing measurements in the outer parts 
to SDSS velocity dispersion determination, using spherical Jeans 
models with a Hernquist profile. All three studies confirmed the 
existence of a mass excess which increases with mass. They all pre- 
ferred a dark matter trend to explain the observations, although they 
could not exclude the IMF variation alternative. 

iGriUo et al.l ( |2009|) compared stellar population masses, de- 
rived from multicolour photometry, to the total masses, inside 
the Einstein radius o f the lenses, published by the SLACS team 
teolton et ai] l2008ah . He assumed for all galaxies the average 
dark matter fra c tion d eterm ined for some of th e galaxies by 
iKoopmans et alj j2006l) and iGavazzi et alj j2007t) . which were 
based on spheric a l mod els a ssuming a s ingle power-law total mass 
and a Hernquist ('1990|) (or|jaffd[l98 ?) stellar profile. From this 
comparison Grillo et al] ( l2009h concluded that elliptical galaxies 
prefer a Salpeter rather than Kroupa/Chabrier IMF. The same result 
was found when comparing stellar dynamical masses for galaxies in 
the Coma cluster by Thomas et al. ( 2009), which includes the dark 
matter halo using fixe d NFW profiles, to stellar population masses 
tarillo & GobaJhoid) . However, given that the dynamical masses 
were determined assuming a fixed NFW profile, the conclusion de- 
pended on the correctness of that assumptio n. A r ecent reanalysis 
of the 16 Coma galaxies by Thomas et al. j201lh . still assuming 
a NFW dark halo profile, interpreted the M/ L excess in massive 
galaxies as more likely due to dark matter, although their results are 
not inco nsistent with a Salpe ter IMF instead. A similar result was 
found bv lWegner et alj j2012l) for 8 galaxies in Abell 262. 

The SLACS team analysed their data assuming all galaxies 
in their sample can be a pproximated by ho mol ogous spher ical and 
is otropic systems w ith a lHernquis^ jl990l) (or |jaffel ll983') profile. 
In lTreu et al] j2010l) they further assumed a spherical NFW profile 
for the dark matter, with fixed slope and only mass as free parame- 
ter. Given that the enclosed total (luminous plus dark) mass inside 
the Einstein radius is essentially fixed by the lens geometry, un- 
der these assumptions the central velocity dispersion is a unique 
function of the profile of the total mass distribution, which in this 
case is defined by a single parameter: the ratio between the stel- 
lar and dark matter components. Comparing the stellar mass from 
the lens mo del to the one from population, based on multicolour 
photometry ('Auger et al."2009^, they fo und that the data prefer a 
Salpeter-like normalization of the IMF. lAuger et al.1 f2010b) used 
the SLACS data and the same spherical Hernquist models but in- 
stead assumed that the halo mass for individual galaxies is known, 
as given by the abundance ma tching techniques, assuming a uni- 
versal IMF (Moster et al."20ltf). They still assumed NFW halos but 
explored both adiabatically contracted ( Gnedin et al. 2004) and not 
contracted halos. They concluded that the data favour a Salpeter- 
like normalization of the IMF over a lighter Kroupa/Chabrier form. 

Although ground-breaking, the analyses of the SLACS data 
depended on some non-obvious assumptions. iGriUo et alj ( l2009h 
had assumed of a fixed dar k matter fraction, as well as a power- 
law form for the total mass. iTreu et alj Jioioh result depended on 
the assumption of a fixed NFW halo slope, as the authors acknowl- 
edged concluding that "the degeneracy between the two [IMF or 
dark matter] interpretations cannot be broken without additional 
information, the data imply that massive early-type galaxies can- 
not have both a universal IMF and universal dark matter halos", in 
agreement with previous dynamical analyses. I Auger et al.1 j2010bl) 
result depended on assuming the halo mass from the abundance 
matching techniques. This led the authors to conclude that "better 
constraints on the star formation efficiency must be obtained from 



the data in order to draw definitive conclusions about the role of a 
mass-dependent IMF relative to CDM halo contraction". 

Other independent studies of the same SLACS data did not 
exclude an IMF variation, but did not co nfirm the nee d for a 
Salpeter IMF to explain the observations. iTortora et alj ( |2010|) 
reached this conclusion u sing spherical Hernquist isotropic mod- 
els like lTreuetal] j2010l) and [AugereTal] j2010bh . but could not 
explain the re ason for the disagree ment. A similar study was also 
performed bv lOeason et all i2012h . They showed that the trend in 
the galaxies enclosed mas ses can be exp lained by a toy model simi- 
lar to the one ad opted by Auger et al (|20_10b), but with a univer- 
sal Rroupa IMF. lOeason etalJ j2012 ) study however did not fit 
the galaxies stellar velocity dis persion, which could be a reason 
for the difference in the results. iBamabe et al. I( l201lh re-analysed 
16 SLACS galaxies using axisymmetric, rather than spherical, dy- 
namical models, and describing the galaxy images in detail, rather 
than assuming fixed Hernquist profiles. Their models are fitted to 
integral-field stellar kinematics, instead of a single velocity disper- 
sion. They find that the data are consistent with bot h a Kroupa or 
a Salp eter IMF, in agreement with the early study bv lFerreras et alj 

boosi) . 

An additional limitation all previous analyses of the full 
SLACS sample was the assumption that all galaxies in their sam- 
ple could be described by homologous spherical Hernquist distribu- 
tions. These models ignore pos sible systematic va r iations of galaxy 
morphology with mass (e.g. ICaon et al] 1 19931 ; iKormendv et alj 
l2009h . Moreover real early-type galaxies, even in the mass range of 
the SLACS survey, are dominated by fast rotating disks (Paper II, 
Paper III), some of which are clearly visible in the SLACS pho- 
tometry. Disk galaxies are not well described by spherical single- 
component models. It is unclear whether galaxy models that do 
not reproduce neither the kinematics nor the photometry of the real 
galaxies under study can be trusted at the ~10% level that is re- 
quired for IMF studies. Possible systematic biases in the results are 
however difficult to estimate. A reanalysis of the excellent SLACS 
dataset using more realistic and flexible models would seem the best 
way to clarify the situation and provide a c onsistent picture. This 
approach has already been demonstrated by Bamabe e t all ( l2012h 
for a single galaxy using the JAM method ( Cap pellari 20081) in com- 
bination with the lensing analysis. One should finally consider that 
lensing studies measure the mass wihin cylinders along the line-of- 
sight. The recovered stellar mass depends on both the assumed dark 
matter profile in the center, as well as at large radii toutton et alj 
l2011iJ) . 

An simple assumptio n on the dark matte r content of early- 
type galaxies was made bv lDuttonetd] ( 12012 b'). They selected the 
most dense galaxies from a large sample of SDSS galaxies and con- 
structed spherical isotropic models to reproduce their stellar veloc- 
ity dispersion. They assumed the total dynamical M/ L of this set 
of galaxies is the same as the stellar one. Comparing the inferred 
mass to the one from stellar population they concluded the galaxies 
require on average a Salpeter IMF normalization. This study still 
depends on the spherical approximation. Moreover it is unclear to 
what accuracy the zero dark matter assumption is verified even for 
the densest galaxies. Nonetheless this study confirms that, unlike 
spiral galaxies, the densest early-type galaxies are not inconsistent 
with a Salpeter IMF. 

In summary, the two result on which nearly all previously dis- 
cussed studies agree are: (i) the total MjL in the central regions 
of galaxies (r < Rc) does not follow the MjL inferred assuming a 
universal IMF; (ii) less massive galaxies, especially spiral ones, re- 
quire and IMF normalization lighter thank Salpeter and consistent 
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with Kroupa-like, while more massive and dense ones allow for a 
Salpeter IMF. Indications were found for the Salpeter IMF to be ac- 
tually required for massive ellipticals dGrillo et al.ll2009l ; lTreu et alj 



l20ld : lAuger et alj2010bl ; lDutton et al.l2012bl) . However these con 
elusions depend on assuming the knowledge of either the dark mat- 
ter fraction or the slope, or the ma ss of t he dark halo. An exception 
is the recent work by iSonnenfeld et al.l j2012l) . which modelled a 
rare elliptical with two concentric Einstein rings, and concluded for 
a Salpeter IMF. More similar objects would be needed to draw solid 
conclusions. Not all studies agreed on the requirement for a Salpeter 
IMF, even when analysing the same SL ACS data ^Ferreras et a l^ 
I2OO8I : iTortora et al]|20ld ISarnabe et al]|201ll : lOeason et allboi C 
It seems that this situation could be resolved by the use of more 
realistic models. 



4.1.5 Systematic IMF variation in galaxies 

A breakthrough in IMF s tu dies was provided by the work 
by Ivan Dokk um & Conrov ( I2OI0I) . further strengthened in 
Ivan Dokkum & Conroy ( 201ll) . Contrary to the previously de- 
scribed set of dynamical and lensing works, they looked for 
evidences of IMF variation in subtle IMF-sensitive spectral fea- 
tures of the near-infrared region of galaxy spectra, in particular the 
Wing-Ford band. In this way their study did not suffer from the de- 
pendency on the halo assumption s. Although the spectral technique 
has been around for decades (e.g. lSpinrad & Tayloilll97l"l) . only the 
availability of reliable stellar population mo dels has made the ap- 
proach sufficiently a c curate for IMF stu dies ( Sc hiavon et al.ll200ol : 
ICenarro et alj|2003l) . Ivan Dokkum & Conrov C201(]h analysed the 
IMF of eight massive ellipticals, from stacked spectra. They used 
new population models that all ow for a variation of the detai led 
abundance patterns of the stars ( Conrov & van Dokkumll2012al) to 
distinguish abundance from IMF variations. They concluded that 
the observed spectra required a bottom-heavy, dwarf-rich IMF. 
Combining their finding with the previous results on the IMF 
of spiral galaxies they tentatively concluded that "Taken at face 
value, our results imply that the form of the IMF is not universal 
but depends on the prevailing physical conditions: Kroupa-like 
in quiet, star-forming disks and dwarf-rich in the progenitors of 
massive elliptical galaxies". 

Our relatively large and well-selected ATLAS^'^ sample and 
high-quality integral-field stellar kinematics appeared well suited 
to resolve the halo degeneracies of previous dynamical and lensing 
analyses and test the claims from the spectral analysis. Contrary 
to previous large studies, we adopted an axisymmetric modelling 
method which describes in detail both the individual ga laxy images 
using the MGE technique CEmsellem et alJll994l : ICapDellari 20oJ) 
and t he richness of our two-dimensional kinematics dCappellaril 
|2008|) . thus avoiding the possible biases of previous more approx- 
imate approaches. For the first time, thanks to the tighter con- 
straints to the models provided by the two-dimensional data, our 
study could leave both the halo slope and its mass as free param- 
eters. The halo slope is allowed to vary in a range which includes 
both the flat inner halos predicted by halo expansion models (e.g. 
iGovematoe t al. 201 0), and the stee p ones predicted by adiabatic 
halo contraction HGnedin et alj|2004l e.g.). Our models also explic- 
itly include the galaxy inclination and anisotropy as free parame- 
ters, although the latter is still assumed to be constant in the region 
where we have data. The parameters are estimated in a Bayesian 
framework with a maximum ignorance (flat) prior on the parame- 
ters. 

We showed that even in this relatively general case the 



models required a va riation of the IMF to reproduce the data 
dCappellari etai]|20I2h . We additionally tested a variety of sensi- 
ble modelling assumptions on the halo, some of which had already 
been employed by previous studies. However this time our mod- 
els accurately described the photometry and kinematics of the real 
galaxies. We confirmed that a non-universal IMF is required under 
all those halo assumption. The similarity is due to the fact that in all 
cases dark matter contributes just about 10-20% to the total mass 
within IRc, so that it cannot explain an observed M/L excess of 
up to a factor 2-3. 

Our ATLAS^'^ sample includes a larger range of masses than 
previous lensing studies. Not only we could quantify the IMF nor- 
malization of early-type galaxies as a class, but we showed that 
the IMF normalization varies systematically within the early-type 
galaxy population, as a function of the stellar (A//L)stars. The IMF 
was found consistent with Kroupa/Chabrier-like at low (i\f/L)stars 
and heavier that Salpeter at large (A//L)stars. This finding bridges 
the gap between the Kroupa-like IMF determinations in spirals and 
the evidence for a Salpeter-like or heavier in early-type galaxies. 
This should be expected given the parallelism and continuity in 
physical parameters between early-type and spiral galaxies as em- 
phasized in our 'comb' morphological diagram (Paper VII). Early- 
type galaxies at low (Af/L)stars in fact completely overlap with 
the region populated by spiral galaxies on the mass-size projection 
of the Virial Plane (Fig. [5). Our result also reconciled the apparent 
disagreement between studies claiming that early-type galaxies in- 
clude cases with Kroupa/Chabrier normalization and those claiming 
that Salpeter is required. 

Although, unsurp risingly, the cleanest IMF trend in 
ICappellari et iD ( |20I2|) was obtained by simply comparing 
the dynamical and population M/L, a trend with velocity dis- 
persion (Te was also presented. This is a natural consequence of 
existence of the M/L — a corr elation (Cappellari et al. 200d; 
Ivan der Marel & van Dokkumll2007,) . In particular Cappellari et alj 
l l20I2h (their fig. 2) showed an IMF trend, with significant scatter, 
going from a Kroupa/Chabrier to a Salpeter IMF within the range 
logj^o(''"e/km s^^) ~ 1.9 — 2.4, with a gradual variation in 
between. 

A flurry of papers have appeared in subsequent months from 
independent groups, all in agree ment on the exi s tence of a sys- 
tematic IM F variation i n galaxies. Spiniello et al.l ( l20I2h used the 
models of IConrov & v an Dokkum (2012a) in combination with 
SDSS spectra and inferred a variation of the IMF from sodium and 
titanium-oxide absorpti ons, which corre late with velocity disper- 
sion (als o illustrated by IZhu et al]|20I(]|) . A similar result was in- 
ferred bv lFerreras et al) ('2012'), also from SDSS spectra, but using 
the population models by Vazdekis et al. (2012). Smith et al. (2012) 
studied the Wing-Ford absorption for galaxies in the Coma Cluster. 
They find that their galaxy spectra are best reproduced by a Salpeter 
IMF and detect a weak IMF trend with metallicity, but do not 
find an IMF trend with velocity dispersion. IConroy & van DokkumI 
( l20I2bl) detected a systematic trend of the IMF, versus either a or 
metallicity, using both optical and near-infrared spectral features. 
A key difference from other spectral studies is that they were able 
to rule out, in a Bayesian framework, a variation of element abun- 
dances as the origin of the empirical trends. 

button et alj ( l20I2ai) showed that it is sufficient to assume the 
knowledge of the halo mass and fit galaxy scaling relations via ap- 
proximate spherical galaxy models, to inf er a systematic variation 
of the IMF. Similarly Tortora et al] ( 1201 2h tested a set of di fferent 
assum ption on the dark halo, using approximate spherical ISersid 
h968l) galaxy models to fit the SDSS velocity dispersion of a large 
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Figure 6. Testing tiie acc uracy of (M/L)salD against tlie completely inde- 
pendent determination of lConrov & van DokkumI j2Q12 b'). which use a dif- 
ferent population code and different spectra with a longer wavelength range 
than the SAURON ones. Assuming all measurements have comparable er- 
rors, the observed scatter between the 35 values is consistent with an en'or 
of 6.5% in (M/L)saip. 



ga laxy sample. Some of t he halo assumptions were also included 
in ICappellari et alj ( l2012h with more detailed dynamical models, 
producing consistent results. Although these results still depend on 
simple approximations and on assumption on the halo, they illus- 
trate that many different halo assumptions, motivated by theoretical 
models, all lead to the same conclusion of a systematic IMF varia- 
tion. 

In summary, after years of debate on whether the IMF or the 
dark matter were responsible for the observed disagreement be- 
tween stellar population and dynamical or tensing indicators, a con- 
sensus seems to have emerged, from dynamical, tensing and spec- 
tral arguments, on a systematic IMF variation with the IMF becom- 
ing heavier with M/L, mass or velocity dispersion. It is time to 
investigate IMF trends against other observables, which we do in 
the following section. 



4.2 IMF variation on the Virial Plane 



In lCappellari et alj ( l2012h we presented a systematic trend between 
the IMF and the stellar {M/L)stais derived using axisymmetric 
JAM models including a dark halo. The IMF was parametrised by 
the IMF mi smatch paramete r a = (M/L)stars/(M/L)saip (in the 
notation of iTreu et alj|201ol) . where (Af/I/)saip was derived from 
full spectral fitting as summarized in Section lZJl Different popula- 
tion approaches were also tested there, producing less clean trends 
but consistent results. The (M/L) saip normalization depends on 
the adopted lower and upper mass cut-offs for the IMF in the popu- 
lation models. The models we use jVazdekis et al.l2oTol) adopt stan- 
dard lower and upper mass cut-ofTs for the IMF of 0. 1 and 100 Mq , 
respectively. (M/L) saip further depends on whether the gas lost 
by the stars during the early stages of their evolution is retained in 
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Figure 7. The systematic variation of the IMF. The (A//L)stars of the 
stellar component, determined via dynamical models, is compared to the 
(M/L)saip determined from spectral fitting using stellar population mod- 
els and assuming for reference a fixed Salpeter IMF. The colours of the sym- 
bols show the galaxies a^, , which was LOESS smoothed in two-dimension 
to emphasize the trends. Diagonal lines illustrate the expected trends if the 
IMF was universal or either the Chabrier (red dash-dotted line), Kroupa 
(green dashed line) or Salpeter (C('ti) m,~^'^, solid magenta line) forms. 
Also shown is the expected trend for two IMFs heavier than Salpeter (blue 
dotted line): a top-heavy one, dominated by stellar remnants (C(™) oc 
m^^'^), and a bottom-heavy one, dominated by dwarfs (Cim) oc m~^ *). 
The black solid curve is a LOESS smoothed version of the data. A clear 
systematic trend is evident, with the IMF being closer to Kroupa/Chabrier 
at the lowest M/L, which also have the lowest ae, and closer to Salpeter or 
heavier at the largest M/ L or largest cTg . A differe nt rendition of this figure 
was presented in fig. 2b of lCappellari et al] j2012l) . 



the central regions we observe or is recycled into stars or expelled 
at large radii. Evidence suggests the gas is not retained in signif- 
icant quantit ies within IJ^e , neither in ionized or hot X-ray emis- 
sion form (Sarzi et alj2010h , nor as a cold component t Young et al] 
hereafter Paper IV). If all the gas was retaine d it would in- 
crease (M/L)saip by up to 30-40% ( lMarastonl200^ . making pop- 
ulation M/L overpredicting total dynamical M/L even for a light 
Kroupa/Chabrier IMF, for a number of galaxies. 

A careful determination of the population M/L for a fixed 
IMF was recently pro v ided fo r 35 galaxies in our sample by 
IConrov & van Dokkum' ('2012b'). They employed different mod- 
els I Con rov & van Dokkum 2012a) from the ones we adopted 
iVazdekis etalj |2010) and used them to fit an independent set of 
spectra, spanning a larger wavelength range than the SAURON ones. 
In Fig.[6]we compare the two determinations of M/L (allowing for 
an arbitrary offset in the absolute normalization, which can differ 
by up to 10% between different authors). Assuming all measure- 
ments have comparable errors, the observed scatter between the 35 
values is consistent with an error of just 6.5% in (A//L)saip. This 
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error i s consistent with the error of 7% predicted bv lGallazzi & Belj 
when optimal spectra information is available. This gives 
confidence that our population (Af/ I/)saip are robust and any trend 
we observe is not due to the details of our population modelling 
approach. 

ICappellari et all JioT^ showed that consistent IMF trends are 
found for a variety of different assumptions on the dark halo, as well 
as for the most general one (with gNFW halo) that leaves both the 
halo slope and mass as free parameters, with only an upper limit on 
the halo slope, deriv ed from model predictions of halo contraction 
dOnedin et alj201 ih (see Paper XIX). The similarity in (A//L)stars 
for the different models is due to the fact that in all cases the data 
allow for a small fraction of dark matter within the region where 
the kinematics are available (^^l.Rc), with the most general model 
implying a median dark matter fractions of just 10%. Given the 
similarity of the different approaches, here we adopt as reference 
the (Af/I/)stars values obtained with a NEW halo, with halo mass 
as free parameter (model (B)). This choice makes it easy to compare 
our results with others that make the same assumption. 

In Fig. [7] we show a di fferent rendition of the similar 
fig. 2b of ICappellari et all ( l2012h . Here we plot (Ai^/I/)saip versus 
{M / L)sta.Ts- We still exclude galaxies with very young stellar pop- 
ulations, selected as having stellar absorption line strength index 
H/3 > 2.3 A. We found that those galaxies have strong gradient in 
the stellar population and this breaks our approximate assumption 
of a constant {M / L)sta,ia within the region where we have kine- 
matics, making both (Af/L)atara and (A^/L)saip inaccurate and 
ill-defined. Different diagonal lines illustrate the expected trends if 
the IMF was universal or either the Chabrier, Kroupa or Salpeter 
{({m) (X m~^'^) forms. Also shown is the expected trend for two 
IMFs heavier than Salpeter: a top-heavy one, dominated by stellar 
remnants (C,{m) oc m~^'^), and a bottom-heavy one, dominated 
by dwarfs (C(?7i) <x m~^'^). This figure clearly illustrates the fact 
that the two indicators of AI/L do not follow a one-to-one rela- 
tion, but deviate systematically, with the IMF being consistent on 
average with Kroupa/Chabrier at the lowest (Af/L)stars and with 
S alpeter or heavier at the largest (Af/L)stars, as already concluded 
in lCappellari et al] ( l2012h . Another way to interpret this plot is by 
noting that galaxies with the largest (Af/L)saip, which are charac- 
terized by the oldest populations, have a Salpeter or heavier IMF, 
while those with the lowest (A//L)saip, which have younger pop- 
ulations, have on average a Kroupa/Chabrier IMF. Detailed trends 
of the IMF with other population indicators will be presented in 
McDermid et al. (in preparation). 

As in the ICappellari et al] ( l2012h we also show with colours 
the galaxies velocity dispersion. However, to emphasize the trend, 
here w e apply the same two-dim ensional LOESS smoothing ap- 
proach dCleveland & Devliij|l988l) introduced in Section[3] instead 
of showing the individual ere values. One can still see the trend for 
the IMF to vary between Kroupa/Chabrier to Salpeter, albeit with 
large scatter, within the interval logj^Q(ac/km s^^) ~ 1.9 — 2.4, 
with a smooth variation in between. This trend will be more pre- 
cisely quantified in Section l431 

In Fig.[8]we show the variation of the IMF mismatch parameter 
(A4'/L)stars/(Af/L)saip On the (A/, ae) and {M,Rc) projections 
of the Virial Plane. The trends in IMF are necessarily noisier that 
the ones in (Af/L)stars, as here the errors in both AI/L combine. 
However the structure in this figure closely resembles the one in 
the previous Fig. [T] |2] and the top panel of Fig. |3] The systematic 
variation in the IMF follows on average the variations in the to- 
tal (Af/L)jAM and its corresponding stellar population indicators 
H/3 and galaxy colour Like the other quantities, also the IMF ap- 



pears to roughly follow lines of constant on the Virial Plane, 
which we showed is tracing the bulge mass at given galaxy size and 
mass. In both projections there is indication for some extra sub- 
structure, with IMF variation at constant tJe . Moreover the proxim- 
ity to the ZOE (thick red line) seems to be an even better indicator 
of IMF than (blue dashed lines). The (Af , R^) projection also 
makes clear why one should expect spirals to have on average light 
IMF: spirals populate the empty zone above our galaxies (Fig. [5] 
and fig. 4 in Paper I) and overlap with the distribution of ETGs with 
the Kroupa/Chabrier IMF. However our plot suggests that, at fixed 
galaxy stellar (or A/jam) mass, the bulge mass (as traced by a^) 
is the main driver of the IMF variation, rather than morphological 
type alone, and one may expect IMF variati ons within the spira l 
population as well. Our finding explains whv lDutton et al.l ( 1201 2bl) 
inferred a Salpeter or heavier IMF when selecting the smallest (or 
densest) galaxies at fixed mass, and assuming no dark matter. Dark 
matter appears to provide a small contribution to the AI/L for our 
entire sample and not just the smallest ones, but the densest galax- 
ies are precisely the ones with the heaviest IMF. A recent claims 
has been made for a Salpeter IMF in the bulge of 5 massive spirals 
jDutton et al.ll20123) . which would go in the suggested direction. 
However the IMF is currently already difficult to infer when the 
Ad/L can be assumed to be spatially constant and the stellar popu- 
lation homogeneous. Relaxing these assumptions makes the results 
more uncertain unless very good population data from galaxy spec- 
tra and accurate non-parametric models for the photometry are em- 
ployed. 

4.3 IMF versus correlation 

In the previous se ction we illustrated the systematic IMF trends al- 
ready presented in lCappellari et al. I( l2012h and we additionally pre- 
sented the variation of the IMF on the projections of the Virial 
Plane. We pointed out that, like the dynamical AI/L, its pop- 
ulation indicators, colour and H/3, and the galaxy concentration 
which traces to the bulge mass, also the IMF broadly follows 
lines of nearly constant tJe on the (A'/jAM,cre) Virial Plane pro- 
jection and consequently constant -y/ GA/jam/(5-R™^'') ~ 
on the (AfjAM, ^e"^'') projection. Interestingly a 'conspiracy' 
between the IMF variation and other galaxies properties is re- 
quired for g alaxy scaling relations to still have a small scatter 
teenzini & C iotti 1993). 

In Fig.[9]we present a direct correlation between the logarithm 
of the IMF mismatch parameter (A//L)stars/(Af/L)saip and the 
logarithm of ae, using the robust linear fitting routine LTSXINEFIT 
presented in (Paper XIX), which explicitly allows and fits for in- 
trinsic scatter and removes outliers. Like earlier we exclude galax- 
ies with very young stellar populations, selected s having stellar 
absorption line strength index H/3 > 2.3 A, leading to a sam- 
ple of 223 out of 260 galaxies. In the fits we quadratically co- 
added JAM modelling errors of 6% (Paper XIX) plus distance errors 
for (A//L)stars, population models errors of 6.5% for (A//L)saip 
(Fig. |6l and 10% for our photometry Paper XXI. We present fits 
for (i) the full sample, (ii) for the subsample of galaxies with SHE 
distances (mostly from Tonry et al. 2001 and Mei et al. 2007: see 
Paper I) and (iii) for galaxies with ae > 90 km s^^, to eliminate a 
few objects with larger scatter. All three relations have a comparable 
observed scatter A « 0.10 dex (26%) in (A^/L)stars/(A//L)saip 
and imply a significant intrinsic scatter of about 20%. When the 
same outliers are removed, we verified that indistinguishable val- 
ues for the parameter, their errors and the in fered intrinsi c scatter 
are obtained with the Bayesian approach of iKellvl ( l2007h . as im- 
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Figure 8. Variation of tiie IMF on tlie Virial Plane projections. Same as in Fig.[T]for the vaiiation of the IMF mismatch parameter (J\//L)stars/(A//L)saip. 
which measure the ratio between the stellar M/ L from dynamical models and the one from population models, with a fixed Salpeter IMF for reference. The 
trends follow lines of roughly constant Ce, or perhaps depends on the distance from the ZOE (thick red line), although some extra structure is visible. 



plemented in his IDL routine LINMIX_ERR. Our three fits provide 
consistent values for the best-fitting slopes, within the errors, and 
nearly consistent normalizations. The fitted relation has the form 



logio 



(M/L)Btars 
(A//L)Salp 



a + & X logj^ 



130km s" 



(6) 



and our preferred values (bottom panel of Fig. |9} have best fit- 
ting parameters and formal errors a = —0.11 ± 0.01 and b — 
0.36 ± 0.06 (parameters and errors for the other fits are given in- 
side the figures). The observed trend of IMF with a appears to ac- 
count for about half of the total trend in the (M/L) — Oe. relation 
(M/L)jAM oc (J°'^^ (see Paper XIX), the remaining one being due 
to stellar population variations. 

Our trend implies a transition of the mean IMF from Kroupa 
to Salpeter in the interval log]^Q((Te/km s~^) ~ 1.9 — 2.4 (or 
(Te ~ 80 — 260 km s~^), with a sm ooth variation in betwe en, con- 
sistently with what can be seen in ICappellari et ^ ( l2012h and in 
Fig. [7] The fact that this trend is slightly weaker than the one im- 
plied by Fig. [7] seems to confirms the intrinsic differences in the 
IMF of individual galaxies. However, part of this difference could 
also be explained if the distance errors were underestimated. One 
way to address this issue would be to repeat our analysis for galax- 
ies in a clusters at intermediate distance like the Coma cluster, for 
which relative distance errors can be neglected. 

Our sl ope is a factor » 3.6 smaller than the "tentative" trend 
reported in iTreu et alj fcoiol) for the same quantities, and making 
the same assumption for the dark halo. A reason for this large dif- 
ference must be due to the fact that their sample only included 
galaxies with a > 200 km s~^ (logjQ((j/km s~^) > 2.3). Our sam- 
ple is too small to reliably study trend in IMF for the galaxies 
above that a^, but their reported trend would exceed the slope of 
the [M/L) — a relation (Paper XIX), making it difficult to ex- 
plain. Their steep inferred trend may then be due to their use of 
spherical and homologous models for all the galaxies, which may 
introduce systematic trends. Moreover they used single stellar pop- 
ulation A'l/L based on colours is expected to be less reliable than 
our determinations based on spectra ( Gallazzi & Bell 2009). A re- 
analysis of the unique SLACS dataset seems required to clarify 
this issue. Our trend with ae is also smaller than the one reported 



by iFerreras et al.l ( 1201 2h from spectral analysis, who find a rapid 
change from a Kroupa to a Salpeter IMF in the narrow interval 
a ~ 150— 200 km s^^. Our relation is not inconsistent with the val- 



ues presented in lConrov & van Dokkuinl(l2012bh also from spectral 
analysis, or with the result reported bv^Spini ello et al.l ( I2012I) . with 
the IMF becoming steeper than Salpeter above a > 200 km s~ ^ . We 
are also broadly consistent with the IMF variation i mplied by the 
non-co ntracted halo spherical dynamical models of iTortora et alj 
l l2012h . Overall there is a qualitative agreement between different 
approaches and the still significant systematic differences in the var- 
ious methods could account for the differences. 



5 DISCUSSION 

5.1 Previous relations as seen on tlie Virial Plane 

We have shown in Paper XVII that galaxy scaling relation can be 
accurately explained by simple virial equilibrium, with galaxies ly- 
ing on a tight Virial Plane (A/jam, o"e, -R™^"), for a large volume- 
limited sample of ETGs (Paper I). Once this has been established, 
the interesting information on galaxy formation is then all contained 
in the distribution and the physical properties of galaxies on this 
plane, which we presented in this paper. 

We find that on the VP: (i) galaxies sizes are delimited by 
a lower-boundary, which has a minimum at a characteristic mass 
MjAM « 3 X 10^° M©; (ii) A number of key galaxy proper- 
ties: dynamical M/L, and its pop ulation indicators H/3 a nd colour, 
which are mainly related to age (Cap pellari et alj|2006h . the nor- 
malization of the IMF and the prominence of the bulge, all tend 
to be constant along lines of constant tTe, on the Virial Plane; 
(iii) Another characteristic mass for galaxy properties is the value 
MjAM ~ 2 X IO^^A/q which separates a region dominated by 
the round or weakly triaxial slow rotators at large masses from one 
dominated dominated by fast-rotators ETGs, flattened in their outer 
parts and with embedded exponential disks (Paper XVII), whose 
characteristics merge smoothly with the ones of spiral galaxies. 
A transition at this galaxy mass appears re quired in our models 
for the formation of fast and slow rotators dKhochfaretallboill , 
hereafter Paper VIII). Interestingly the approximate stellar mass 
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Figure 9. IMF versus (Je coirelations. The plots show as blue filled 
circles with error bars the logarithm of the IMF mismatch parameter 
(M/L)stars/(Af/L)saip versus the effective velocity dispersion a^, for 
the subsets of 223 ATLAS^° galaxies with H/3 > 2.3 A. Green symbols 
are outliers automatically removed from the fits. Dashed and dotted red lines 
indicate the Icr and 2.6cr (99%) observed scatter around the best-fitting re- 
lation (black solid line). The top panel includes all 223 galaxies. The middle 
panel shows galaxies with accurate SBF distances. The bottom panel shows 
galaxies with > 90 km s~^. 



« 1 X 10 Mq where normal early-type galaxies starts to 
appear corresponds to the threshold for q uenching of field galaxies 
recently discovered bv lGehaetal1 ( l2012i) . Below that mass only the 
cluster or group environment can strip spirals of their gas. 

Although our sample is limited to a minimum mass 
MjAM^O X 10^ Mq, our picture extends smoothly to lower 
masses. A comparison of the trends outlined in fig. 4 of Pa- 
per I or Fig. |5] here, with si milar rela tions for lower mass galax- 
ies (e.g. fig. 7 of Binggeli et al | [ l984l : fi g. 38 of Kormendy et all 
l2009l : fig. 12 of'chen et al. 2010; fig. 4 of Misgel d & Hilkeill201 ll 
fig. 20 of lKormend v & Bender 2012) clearly shows that our ETGs 
trends continues with the dwarf spheroidal (Sph) sequence at lower 
masses, while the spirals sequence continues with a sequence 
of low-mass l ate spirals or irregula r s (Sc- Irr), as independently 
noted also by iKormendv & Bendeij ( l2012h . The parallelism be- 
tween ETGs and spiral galaxies in scaling relations was one of the 
driver for our proposed revision (Paper VII) to the tuning-fork di- 
agram bv lHubblel l ll936l) , in the spirit o f Ivan den Bergl] Jl97d ). In 
Fig. [To] we combine our results on morphology, kinematics and 
scaling relations in a single diagram, using the same galaxy symbols 
as Paper VII. This picture allows us to provide a new perspective 
and a clean empirical view of a number of classic scaling relation 
and known trends in galaxy properties. 

The classic Faber & Jackson relation between a — L 

is well known to be a projection of the FP. We study it here using 
mass instead of light, and find that it is tracing the envelope de- 
fined by the ZOE. We find that the relation is not well described by 
an approximately linear trend. This is consistent with some earlier 
claims of a possible change in the slope o f the (t—L relation of ellip- 
tical galaxies at the low l uminosity end ( lTonrvl[l"98ll : lDavies et alj 
|1983| J Tortora et ai] |2009h . However even recent observations from 
the large SDSS sample have subsequently failed to find evidence 
for a cle ar curvature, likely due to the insufficient quality of the data 
teemardi et al. 2003 ; Gallazzi et al..2006.) . For this reason the rela- 
tion is generally assumed to be a power-law, except in the region of 
dwarf elliptic als (de Riicke et al."2 005l ; lMatkovic & Guzman(l2005l ; 
iForbes et al.f 2008; Cody et al. 200^. The reason for this confusion 
is that the (A/jam , projection is not too far from an edge-on view 
so that the bend in the ZOE is not well pronounced and the evidence 
for the bend can be robustly detected here for the first time in ETGs. 
Samples of morphologically selected elliptical galaxies, for which 
the relation is usually derived, tend to populate mainly the top (red) 
region of the diagram, where the effect of the cusp is not evident. 
Moreover most previous studies did not reliably sample the lowest 
a regime. 

Another well known projection of the FP is the iKormendvl 

relation. When using mas s instead of light i t beco mes clear 
it represents the analogue of the iFaber & Jacksonl l ll976l) . but this 
time in the (A^jam, -R™^'') projection of the VP. Also in this case, 
when samples are morphologically selected to consist of ellipticals, 
they tend to populate mostly the region of the diagram near the 
ZOE, defining a relatively narrow sequenc e (IKormendv et al .l2009i ; 
IChen et aljlioicl; iMisgeld & IIilkeJl201 ih . Although the sequence 
is useful for a number of studies, it is important to realize that it is 
not a real sequence in galaxy space on the VP. It is due to the sample 
selection and it represents essentially one of the contour levels of a 
continuous trend of galaxy properties, spanning from spiral galax- 
ies, to ETGs, and only terminating on the well defined ZOE (Fig. [5] 
and figure 4 in Paper I). 

The characteristic mass MjAM ~ 3 X l() ^ °Afa is the same 
transition mass discovered by iKauffmann et al.l (l2p03a) who state 
that "low-redshift galaxies divide into two distinct families at a stel- 
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Figure 10. Schematic summai'y of the resuhs presented in Section[3]and Section|4] ETGs properties, dynamical M/L (Fig.[T) or its population proxies, H/3 
and galaxy colour (Fig.|2), concentration (Fig.|3), which traces bulge mass, and IMF mass normalization (Fig.[8), all tend to vary along lines of nearly constant 
yjGhUK^R^^^) o-e. This sequence of ETGs properties merges smoothly with the one of spiral galaxies, with little overlap between late spirals (Sc-Irr) 
and ETGs, a significant overlap between early spirals (Sa-Sb) and fast-rotator ETGs with low M/L and no overlaps between spirals and fast-rotators with high 
M/L. Three characteristic masses are emphasized in this diagram: (i) below Ri 10^ Mq there are no regular ETGs and the mass-size lower boundary 
is increasing; (ii) M* Ri 3 X 10^0 Mq is the mass at which ETGs reach their minimum size (or maximum stellar density), before a sudden change in slope 
Rc Qc M^-"^^ at larger masses (see also fig. 4 in Paper I); (iii) Below M-^ 2 X 10^^ Mq ETGs are dominated by flat fast rotators, showing evidence for 
disks (Paper XVII), while slow rotators are rare. Above this mass the population is dominated by round or weakly triaxial slow rotators. These smooth trends 
in scaling relations motivated our proposed parallelism between spirals and ETGs. To emphasize this connection we uses the same morphology symbols as our 
'comb' diagram in Paper VII. 



lar mass of 3 X 10^° Mq. Lower -mass galaxies have young stellar 
populations, low surface mass densities and the low concentrations 
typical of discs. Their star formation histories are more strongly cor- 
related with surface mass density than with stellar mass." A simi- 
lar trend involving colours and also better correlated with surface 
density (or with the velocity dispersion inferred from the scalar 
virial relation) than with mass, was found to extend to redshift up 
to 2 ~ 3, with red galaxies being s ystematically smal l, and blue 
galaxies being large at a given mass jFranx et alJboOSh . This was 
recent ly confirmed, still using photometric data alone, bv lSell et alJ 
The correctness of all these statements can now be easily 
and accurately verified for the ETGs subset in Fig.[T}{3] 

The novelty of our work, with respect to all previous stud- 
ies, is that we have unprecedentedly accurate and unbiased dy- 
namical masses, and stellar velocity dispersions, instead of inferred 
value. This allow us to conclusively state that neither dynamical 
mass nor stellar surface density are actually the best descriptor 
of galaxy properties, the main trend being along the cTe direction 
(which includes rotation in the case of di sk galaxies ) . Our c lean 
ATLAS^° result was already presented in iCappellarilfcoilah and 
subseque ntly confirmed wit h SDSS data and using virial mass esti- 
mates bv lWakeetal] ( l2012l) . The trend we observed for the ETGs 
can be extended to spiral galaxies, which fill the region of larger 
sizes above the ETGs in the (A/jam, -R™^'') projection, smoothly 
overlapping with the ETGs for the largest spiral bulge fractions (fig- 
ure 4 of Paper I and Fig. [5](- It has been known for long that in 
spirals luminosity-weighted ages are lower, star formation is larger 
and colo urs bluer on average than the ETGs , essentially by defini- 
tion (e.g. lHubblelll93d Ivan den Berghll 19761) . 

Similar trends between age and size (or surface brightness), 
with older o bjects being smaller at giv en age, were found to per- 
sist in ETGs. Ivan derWeletalJ ( l2009al) state that "at a given stellar 
velocity dispersion, SDSS data show that there is no relation be- 



tween size and age". The s ame age-size trend was pointed out by 
.Shankar & Bemardil ( I2OO9I) . and in different terms bv Graves et alj 
l l2009l) who state that "no stellar population property shows any 
dependence on Rc at fixed cr, suggesting that a and not dy- 
namical mass is the better predictor of past SFH". The age-size 
trend was also c o nfirm ed in different samples of g alaxies by 
IValentinuzzi et alj ( |201(]|) and iNapolitano et al.l ( l201(j) . All these 
findings are another way of saying that age variations must fol- 
low lines of constant a on the VP as we find here for M/L, Hj3 
and colour, and confirm using age for o ur sample in McDerm id et 
al. (in preparation) (see also figure 15 of Gallazzi et al. 2006). The 
only contrasting view is the one bv iTruiillo et al.. (2011.) , who find 
a lack of age-size trend both at low and high-redshift. 

The same characteristics mass MjAM « 3 X Mq of 

iKauffmann et al.l (^0Q3a), which constitute the locatio n of the break 
in our ZOE, was found bv lHvde & Bemardil ( [2009al) in the mass- 
size relation of 5 x 10* SDSS galaxies, their trend is significant 
but quite subtle. The reason for the curvature they find becomes 
clear from what we find: the average radius of the objects on the 
VP at constant mass, bends upwards at low masses, due to the cusp 
in the ZOE, with the strength of the effect being dependent on the 
specific criterion adopted to select ETGs. The bend in the mass- 
a however is essentially undetected in the SDSS data, likely due 



to the insufficient quality 
l lHvde & Bemardill2009 j) . 



A curvature i n the l uminosity- size relation was presented in 
iGraham & WorlevI ( l2008t) (see also Ichen et"ai]|2010l) . It appears 
when joining dwarf spheroidal galaxies (Sph) and el l iptical s as 
a continuous sequence, with smoothly varying ISersid j 19681) in- 
dex. The observed continuity of the relation is interpreted as 
a physical connection be tween dwarf sphero i dals a nd ellipticals. 
The tr end discussed by Graham & Worlev 1 20()8|) can also be 
seen m iBinggeU et al] ( Il984h , iKormendvl "( ll985l) . lKormendvetalJ 



© 2012 RAS, MNRAS 000,[Tll25l 



The ATLAS^^ project - XX. Mass-size 1 9 



jlOOSl) and lMisgeld & HilkedlOlll . where the dwarf spheroidals se- 
quence appears to sharply bend from the elhpticals sequence. But 
iKormendvl ( IT983) interprets dwarf spheroidal as constituti ng a sep- 
arate fami ly, of gas-stripped dwarf sp irals/irregulars (see iGrahamI 
I2OIII and iKormendv & Bended 1201 2| for two complementary re- 
views of this subject). 

Our results cannot not be compared to theirs in a statisti- 
cal sense, as galaxies in their diagrams are, by design, not rep- 
resentative of the population in the nearby Universe, and cer- 
tain classes of objects (e.g. M32) are by design over-represented. 
Our sample is volume-limited and for this reason gives a statis- 
tically representative view of the galaxy population above a cer- 
tain mass. Still the fact that the sequence of dwarf spheroidals 
and low-mass spirals/irregulars, lie on the continuation of our 
trends for fast rotator ETGs and spiral galaxies respectively, be- 
low the M<6 X l(f Mq mass limit of our survey, suggests 
a continuity between dwarf spheroidals and the low-mass end 
of our disk-dominated fast-rotator ETGs population, which in 
turns we showed are closely related to spiral galaxies. For this 
reason our results reconciles the apparent contr a st between the 
findings of a Sph-E dichotomy tKormendv 19851: Ide Riicke et 



I2OO5I : Ijanz & Lisked I2OO8I: iKormendv et al 



20091) and the ones 



of a continuity taraham & GuzmaiJ I2063F loavazzi et alj l2005l : 
iGraham & Worie v 2008; Forbes et al.) |201 ih . Our finding in fact 
agrees with the proposed common origin of dwarf spheroidal 
and low-mass spir al galaxies and irregulars jKormendvl 1 19851 : 
lDekel&Silklll986h . but also shows an empirical continuity be- 
tween dwarf spheroidals and a subset of the ellipticals family. The 
missing link betw een Sph and E is constituted by disk-dominated 
fast rotator ETGs ( lEmsellem et alJl2007l : ICapDellari et alj|2007t) . In 
fact the continuity is not between "true" ellipticals, namely the 
slow rotator, and dwarf spheroidals, but between "misclassified" 
ellipticals with disks and SO, namely the fast rotators, and dwarf 
spheroidals. After this text was written and the parallelism be- 
tween spiral galaxies and early-type were presented (Paper VII) to 
interpret t he trends we obse rved in galaxy scaling relations (Pa- 
per I and ICappellaril 1201 lah . a confirmation of this picture, in- 
cluding its extension to low mass systems was also provided by 
IKormendv & Bender (2012J. 

In summary our view of the VP is consistent and provides a 
confirmation and the cleanest empirical view of all these reported 
trends between galaxy properties, which we interpreted as due to 
the variation of the bulge fraction in fig. 4 of Paper I and explained 
as due to an apparent conversion of spirals into fast rotators ETGs 
in Paper VII. 



5.2 Implications for galaxy formation 

Galaxy formation is the superposition of a number of complex 
events that happen in parallel. Here we sketch a tentative picture 
of some of the phenomena t hat can play a rol e in explaining what 
we see. We refer the reader to lBois et alj ( l201 ll hereafter Paper VI), 
Paper VIII and Naab et al. (in preparation) for a more in-depth dis- 
cussion. 

The smoothness and regularity of the trends we observe and 
the fact that they extend to spirals, appears to indicate a close con- 
nection between the formation of the two classes of objects (Pa- 
per VIII). The same similarity between the fast rotator ETGs and 
spirals in terms of their morphology and degree of rotation, lead us 
to propose a revision (Pa per VII) of the classic morphological clas- 
sification ( lHubble|[l93db to emphasize the p arallelism between the 
fast rotators and spirals, in the same way that lvan den BerglJ ( ll976h 



proposed it for SOs and spirals. Only the most massive slow ro- 
tators appear to form an empirically separated class. The kinematic 
morphology-density relation (Paper VII), which apphes to our kine- 
matic classes the relation discovered by iDressled l [l980b . suggest 
that most spirals are being transformed into fast rotators due to en- 
vironmental effects (Khochfar & Ostriker 2008), with a mechanism 
that is sufficiently 'gentle' to preserve the near axisymmetry of the 
disk (Paper 11). 

One process which is often mentioned in the context of clus- 
ter galaxy populations is stripping by the interstellar medium 
I S pitzer & B aade 1951; Abadi et al. 1999), which may impact the 
global morphology of the galaxy by removing a significant fraction 
(or most) of its gas content but mostly preserving the stellar disk 
component. This may explain some of the most flattened galaxies in 
the fast rotator class (e.g., NGC4762), as emphasised in Paper VII. 
This would work by removing some of the mass of the galaxy keep- 
ing a relatively constant effective size (of the stellar component), 
and may contribute to the scatter in the Mass-size plane for fast ro- 
tator as well as to their overlap with the more gas-rich spiral galax- 
ies. 

As to explain the intermediate to high mass end of the fast 
rotator class, we would need a process which is able to increase 
the bulge size, while at the same time removing the gas or shut- 
ting off star formation (Paper VIII). The empirical signatures of 
this phenomenon are visible in our data as an 'apparent' decrease 
of the galaxy size, and increase of CTj, which is actually due to the 
bulge concentrating more light at smaller radii, accompanied by a 
increase in M/L, which is tracing a decrease of star formation or an 
age increase. The process appears to generally preserve the intrinsic 
flatness of the stellar disks at large radii (mid dle panel of Fig .[3]l. 
A similar scenario was recently proposed by iBell et all ( I20I2I) to 
interpret the relationship between rest-frame optical colour, stellar 
mass, star formation activity, and galaxy structure from 2; ~ 2 to 
the present day. 

Intense gas-r ich accretion events, mostly via cold streams 
jPekel etal.ll2009l) . or major gas rich mergers (Paper VIII), will 
increase both the mass an d Sg. During the accret ion the gas 
may sink toward the centre jMihos & Hernquisl|[l994l) until it be- 
comes self-gravitating and starts forming stars. It is during this 
phase of rapid gas accretion that the M/L ~ a jCappellari et alj 
'200^ and IMF - cr (Section gjU relations, th e tilt of the F P 
(Dressier et al. 1987; Faber et al. 1987; Diorgovski ~& Davislll987h . 
theMgfc-o- (Burstein et al. 1988; Bender et al. 1993) or the Mg6- 
T4sc relation dDavies et al., ,1993. ; .Scott et aU ,20091. Paper XXI) 
will be imprinted in the ETGs population 1 Robertson et alj|2006l ; 
lHopkinsetai]|2009bh and then mostly preserved in the following 
evolution. 

The early progenitors of today's fast rotator ETGs would 
be high-redshift spirals, which are different from local one s. 



They are characterized by giant gas clumps (Elmegre en et al 
'2007; Genzel et al. 2011') have high gas fractions (T acconi et al 



i2010; Daddi et al.. .2010) and possess large veloc ity dispersion 
jFQrster Schreiber et alJl2006ll2009l ; rLaw et alj20il) . In that situa- 
tion bulges may form naturally as the clumps collide and may either 
sink to the centre (Bournaud et al. 2007) or are efficiently de stroyed 
by stellar feedback |Gen el et al. 2012; Hopkins et al.ll20ll]) . Secu- 
lar effects jKormendv & Kennicutt.20o4) will also contribute to the 
bulge growth and Ee increase, while keeping the mass unchanged, 
as will contribute the destabilizing effect of minor mergers. 

During the bulge growth some process must be able to turn 
off star formation, without destroying the fast rotating disks that 
still dominate the local ETGs population (Paper II; Paper III) 
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Figure 11. Evolution scenario for ETGs. Tlie symbols are the same as in 
Fig. [T] while the large aiTows indicate the proposed interpretation of the 
observed distribution as due to a combination of two processes (a) bulge 
or spheroid growth, which seems associated to the quenching of star for- 
mation, which moves galaxies to the right of towards the bottom, due to 
the increased concentration (decreasing iJc), and (b) dry merging, increas- 
ing Rc by moving galaxies along lines of roughly constant Ue (or steeper), 
while leaving the population unchanged. A schematic illustration of these 
two processes is shown in Fig. 1121 



and that dominates the ETGs pop ulation already from z ^ 2 
jvan der Wei et al.l 1201 ll : see also Ivan Dokkuri] llOllh . The re- 
duced efficiency o f star formation by morphological quenching 
dMartig et ai]|2009h may be one of the processes explaining why 
on average bigger bulges correspond to older ages and larger Af/L. 
But bulge/spheroid growth seems also associated with AGN feed- 
back, which would provid e another bulge-related quenching mech- 
anism jSilk&Reesl ll998'). During the sequence of bulge/spheroid- 
growth followed by quenching, the original gas rich spiral will 
move from the left of the observed (Mjam, i?™'"') plane towards 
the right, or possibly the bottom right due to the increased con- 
centration, while intersecting the constant a lines (blue arrow in 
Fig. II It. At the end of the evens the galaxy will be a fast rotator 
ETG, generally more massive and with a bigger bulge (smaller i?c, 
larger and concentration) than the precursor clumpy spiral. 

At any stage during the bulge-growth phase the galaxy may ac- 
crete purely stellar satellites iKhochfar & Silkl2006h . In the case of 
these dry mergers the situation is quite different and one can predict 



the fi n al configuration using energy conservation jHemauist et al. 
19931: ICole et alj I2OO0I : iBovlan-Kolchin et al.l I2OO6I : ICiotti et al. 



20071) . The predictions show that sizes will increase as the mass 



grows. For major mergers (equal mass) and typical orbital config- 
urations, one can show that the mass and radius dou ble, while a 
rem ains constant in ag ree ment with simulati ons CNipoti et alj2009l) 
(see iHilz et aLl l2012bl and iHilz et al J 120123 for more accurate nu- 
merical simulations and detailed physical explanations of this pro- 
cess). While in the limit in which the same mass doubling hap- 
pens via small sa tellites, as mostly expected from the shape of the 
ISchechteil ( fT973) mass function, the radius will increase by a factor 



Figure 12. Schematic representation of the two main processes responsible 
for the formation of the observed distribution of galaxies on the VP. (a) bulge 
growth via cold accretion, secular evolution, or minor mergers, followed by 
quenching by AGN or other mechanisms, leaving the galaxy more massive, 
more compact, and consequently with a larger u^, and gas poor (blue arrow 
in Fig. llU ; (b) major or minor dry mergers, increasing galaxy mass and sizes 
at nearly constant o-g, or with a possible decrease, leaving the population 
mostly unchanged (red arrow in 

Fig.QD- (taken trom .Cappellari 201 Ibh . 



of four and the dispersion will be twice smaller 1 Naab et al. |2009| : 



iBezanson et al.l2009l : lHopkins et alj2009j : lHilz et alj2012bl) . As a 

result the galaxy will move along lines that are parallel to the con- 
stant a lines, or steeper. During these dry mergers, given that there 
is little gas involved, the stellar population, colour and MjL will 
remain unchanged (red arrow in Fig.IIB. 

The gas-rich mergers/accretion scenario is generally consistent 
with the observed correlation between supermas sive black holes 
(BH ) and galaxy velocity dispersion or bulge mass jPi Matteo et alj 
|2005|) . It is generally believed that the correlations indicate a 
joint evolution of galaxies and BHs, with BH growth happen- 
ing at the same time as the bulge growth , and proyiding_a self- 
regulation via feedback jSilk & Reel 19981 : but see |Pendl200 7' and 
Jahnke & Maccio 201 1 for a non-causal origin of the correlations). 
We demonstrate that in deed the a vari ation directly traces the bulge 
growth iKormendv & Gebhardlll200l1) . In fact a traces the central 
galaxy concentration and the bulge size as estimated on optical mor- 
phology, while the outer galaxy disks remain flat. This implies that, 
if the BH indeed accretes from the same gas that grows the bulge, 
BH ma ss should correlate better with a than with total mass as ob- 
served taebhardt et alJl2000l : lFerrarese & Menit^lioOQl 

However, when ETGs experience dry mergers, their BH grows 
in proportion to the mass, but galaxies move along lines of nearly 
constant a. For this reason, an expectation from this picture is 
that BHs at the high mass end should start to appear too massive 
with respect to the predictions of the Mbh — o relation (see also 
iNipoti et alJl2003l : lBovIan-Kolchin et ai]|2006l) . The high-mass BH 
end is still not sufficiently populated to reliably test this prediction, 
but in direct evidence seems to support this possibility (.Lauer et al] 
l2007l) . as well as the recent detection of two giant black holes at 
the centre of two bright cluster galaxies, which are clearly more 
massive than the A/bh — o prediction (McConnell et al. 201 1). 

In summary we propose that the distribution of galaxy prop- 
erties on the VP, where MjL and age follows lines of constant 
a on the VP, could be explained by the combination of two pro- 
cesses, which can happen multiple times during the evolution of a 
single galaxy: (i) They accrete gas, which grows the bulge and BH, 
shrinks their i?c, and increases and concentration, while some 
process which seems associated with the bulge growth (e.g. AGN 
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feedback), quenches star formation; (ii) They experience mostly mi- 
nor dry mergers that move them along lines of roughly constant o 
(or steeper). 

An open question in the scenario in which ETGs evolve rela- 
tively quietly from spirals comes from the comparison of our find- 
ings with the empirical scaling relations one observes at larger red- 
shift. In fact at z > 1.5 gal axies are found to be smaller than local 
ones with the same mass toaddi et all 12001 i Truiillo et alj|2006l . 
I2OO7I : Ivan Dokkum et al.l I2OO8: Cim atti et alj boOSh . They popu- 
late the region below the ZOE of local galaxies. This may indi- 
cate that the compact high-redshift ETGs follow a different and 
more viole nt evolutionary pa th than the more quiet majority of lo- 
cal ETGs tearro et al.ll2012b . We show i n the to p panel of Fig. [2] 
that the 12 available de te rminations of a jCenarro & Trujillo 2009; 
Cappellari et al.' '2009'; 'van Do kkum et alj 120091; lOnodera et alj 
2010 ; van de Sande et al. 201 1) are all consistent within the er- 
rors with the ATL AS^° VP, with the only e xception being the 
object presented in Ivan Dokkum et all ( l2009h . A way to recon- 
cile this very mild (or lack of) a evolution is by assuming that 
the compact primordial ETGs grow mostly by accretion of small 
satel lites in their outer halos, while preserving the central struc - 
ture jNaab et alj200^ ; lHopkins et al.l2010l ; IOser et al.l201oll2012l) . 
This seems consiste nt with the shape of the ph otometric profiles 
of the early ETGs (iHqpkin s et al. 2009a; Bez anson etaU |2009| ; 
Ivan Dokkum et al.ll2010l ; iHilz et alj|2012al) . A caveat is that sig- 
nificant biases may still exist in the high-redshift photometry 
jMancini et allbOlOl) . considering that systematic differences of up 
to a factor of two exists even on well observed ETGs in the nearby 
Universe (Kormendv et al. 2009, Chen et al. 2010, Paper I). More- 
over comparisons of photometric profiles tend to be made against 
bona fide ellipticals, while the remnant of the high-redshift ETGs 
are likely disk-like fast rotators ETG, which have systematically 
different profiles. Finally the comparisons should ideally be done 
in mass density, instead of surface density, but kinematic informa- 
tion is available for only a handful of galaxies. This implies that 
there is perhaps still some room for the compact high- 2: ETGs to 
become more consistent with local ones, than currently assumed. 



6 SUMMARY 

In the companion Paper XIX we describe in detail the dynami- 
cal models for all the ATLAS ETGs that were introduced in 
ICappellariet"ai] l l2012h . We found that galaxies lie on a thin Virial 
Plane (VP) in the three-dimensional parameter space defined by 

(MjAM,ae,i?r")- 

Here we studied the inclined projection of the VP and find 
that: (i) the location galaxies define a clear boundary, well de- 
scribed by two power-laws, joined by a break at a characteristic 
mass MjAM = L x {M/L)e « 3 x IO^^Mq, which corresponds 
to the regime where a number of global galaxy properties were re- 
ported to change; (ii) The distribution of {AI/L)e, as well as pop- 
ulation indicators of [M / L)^op like HjS and colour, together with 
galaxy concentration that we show here is measuring bulge size, all 
tend to be constant along lines of constant CTe. 

These findings explains why CTe (not L, M, RcOT Ee) has of- 
ten been the most successful single descriptor of galaxy properties. 
We discuss a number of previously found galaxy correlation and we 
find that they are all consistent with what we find. In fact most ob- 
served relations turn out to be just special projections of the cleaner 
and more general view provided by the VP. 

The stellar initial mass function was previously shown (e.g. 



ICappellari et alj|2012l) to vary systematically with (i\f/L)stars and 
as expected it also follows the other population indicators and also 
tend to vary along lines of nearly constant on the VP. This is a 
necessary expectation if one wants to preserve the tightness of the 
Fundamental Plane or the (M/L) ~ tJe relation. The trend of IMF 
with a appears to account for about half of the total trend in the 
(M/L) — (Te relation (see Paper XIX), the remaining one being due 
to stellar population variations. 

The distribution of galaxy properties on the VP can be qualita- 
tively interpreted as due to a combination of two main processes: (i) 
bulge growth, changing the galaxy population and decreasing _Rc, 
while increasing cTe, and (ii) dry merging, increasing by moving 
galaxies along lines of approximately constant a^, while leaving 
the population unchanged. It is unclear where the reported dense 
ETGs fit in this picture. They may simply follow a different route, 
which forms the most massive galaxies and is not affecting the fast 
rotator population, which dominates ETGs in the nearby Universe. 
Or there may be some remaining biases in the photometric observa- 
tions, producing an underestimation of the radius. The few available 
high-redshift a determinations are in fact all with one exception 
consistent with what we observed on the VP. Stellar kinematics of 
high-redshift ETGs is critically needed to address this question. 
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